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FOREWORD 


This  report,  describing  a  fluid-mechanical  model  of  the  thoraco¬ 
abdominal  system  derived  using  theoretical  concepts  and  empirical 
data,  incorporates  dimensional  analysis,  physical  and  biological 
scaling,  and  evaluation  of  critical  parameters  to  help  elucidate  mammal¬ 
ian  response  to  sudden  changes  in  environmental  pressure. 

Though  the  work  shares  many  limitations  common  to  "state  of 
the  art"  theoretical  explorations,  it  is  nonetheless  pertinent  to  the 
important  area  of  estimating  human  responses  to  blast  overpressures 
from  those  measured  and  known  to  be  hazardous  in  several  species 
of  animals  of  different  sizes. 

The  theoretical,  computational  and  empirical  data  included  in 
the  present  study  represent  part  of  an  ongoing,  long-term  program 
aimed  at  clarifying  the  biological  effects  of  blast-induced  phenomena 
and  assessing  the  consequences  of  exposure  thereto. 
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ABSTRACT 


A  mathematical  model  was  described  which  was  developed  to 
compute  some  of  the  fluid-mechanical  responses  of  the  thoraco¬ 
abdominal  system  subjected  to  rapid  changes  in  environmental 
pressure'.  The  approach  --  helpful  in  understanding  many  of  the 
primary  effects  of  air  blast  on  animals,  but  applicable  to  related 
problems  involving  rapid  changes  in  environmental  pressure  as 
well  --  incorporated  an  air  cavity  representing  the  gaseous 
volume  of  the  lungs  (although  it  is  recognized  that  gas  in  the 
organs  of  the  abdomen  may  influence  the  response  of  the  system), 
two  movable  pistons  and  an  orifice  through  which  gas  might  pass 
in  either  direction.  One  of  the  pistons  represented  the  chest  wall 
and  the  other  that  portion  of  the  abdomen  which  moves  with  the 
diaphragm  to  change  the  lung  volume.  Each  piston  was  "assigned" 
an  effective  mass  and  area,  a  spring  constant,  and  a  damping  factor. 
The  orifice  was  taken  to  "incorporate"  the  characteristics  of  the 
many  airways  of  the  respiratory  system.  Parameters  relating 
the  animal  to  the  model  were  estimated,  tested  and  then  adjusted 
as  required  by  comparing  model  results  with  experimental  records 
of  thoracic  pressures  recorded  for  rabbits  exposed  to  blast  waves 
in  shock  tubes. 

Equations  were  derived  to  scale  parameters  applicable  to  a 
given  animal  to  those  for  similar  creatures  of  arbitrary  mass. 

I  f  dimensional  analysis  other  equations  were  developed  to  relate, 
for  a  given  biological  response,  the  body  mass  of  similar  animals 
to  blast  wave  parameters.  Numerical  solutions  of  the  model  were 
presented  to  help  explain  the  mechanisms  involved  when  animals 
were  "loaded"  \&ith  typical  wave  forms  or  with  pulses  increasing 
to  a  maximum  in  a  stepwise  manner,  a  contingency  associated 
with  a  quite  significant  increase  in  mammalian  tolerance  to  over¬ 
pressure.  Differences  in  response  to  "short-"  and  "long" -duration 
blast  waves  were  noted.  Applications  of  the  scaling  concepts  were 
exemplified  in  several  ways  making  use  of  the  published  data  in 
blast  biology.  In  one  instance,  the  blast  tolerance  of  a  70-kg 
mammal  was  estimated  for  sea-level  ambient  pressure  making  use 
of  experimental  data  for  dogs  and  goats  obtained  at  Albuquerque  al¬ 
titude  (ambient  pressure  of  12  psi).  That  the  latter  might  have 
significant  implications  in  assessing  human  response  to  blast- 
produced  overpressures  was  discussed  along  with  several  other 
relevant  matters. 
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A  FLUID-MECHANICAL  MODEL  OF  THE  THORA  CO -ABDOMINAL 
SYSTEM  WITH  APPLICATIONS  TO  BLAST  BIOLOGY 


1.0  INTRODUCTION 

In  the  course  of  studies  in  blast  and  shock  biology  over  the 
past  decade,  the  investigations  concerned  with  the  consequences  of 
exposure  to  blast-induced  variations  in  environmental  pressure 
(primary  blast  effects)  have  been  given  increased  attention  for  at  least 
three  reasons.  The  first  is  related  to  the  fact  that  early  lethality  was 
known,  and  further  demonstrated,  to  be  among  the  characteristic  re¬ 
sponses  for  all  species  of  mammals  exposed  to  appropriate  rises  in  air 
pressure  produced  by  actual  detonations  or  by  simulated  explosions, 
i.  e.  ,  by  means  of  shock  tubes.  The  second  reason  is  that  the 

magnitude  of  a  hazardous  overpressure  is  significantly  lower  for 
"long" -duration,  "fast"-rising  wave  forms  that  can  be  produced  by 
nuclear  explosives  than  it  is  for  similar  blast  waves  of  shorter 
duration  resulting  from  detonations  of  conventional  charges  of  practical 
sizes. •**  6*  7-1  IT  The  third  reason  is  related  to  the  experimental 
observation  that  mammals  are  extraordinarily  sensitive  to  the  rate  and 
character  of  the  pressure  rise.  *0-13,  15,  18-21  por  eXample, 
an  otherwise  lethal  blast  wave  car  _>e  tolerated  if  there  is  a  sufficient 
delay  in  the  time  to  maximal  pressure  or  if  the  pressure  is  applied 
in  two  fast-rising  steps  separated  by  an  almost  unbelievably  short 
interval  in  time;  viz. ,  several  tens  to  hundreds  of  microseconds 
depending  upon  animal  size.  10*  18-23 


Though  damage  to  the  lung  and  the  appearance  of  gas  ' 
emboli  in  the  coronary  and  other  arteries  in  blasted  mammals  no 
doubt  help  explain  the  characteristic  early  lethality  mentioned 
above,  3~->>  7-10,  1 5-19,  22-26  the  detailed  events  which  are  respon¬ 
sible  for  the  origin  of  the  hemorrhages  and  emboli,  until  recently, 
remained  obscure  as  did  a  satisfactory  understanding  of  the  variations 
in  tolerance  alluded  to  previously  as  connected  with  changes  in  the 
pulse  duration  and  the  rate  and  character  of  the  pressure  rise.  These 
matters  emphasized  the  need  for  paying  careful  attention  to  all  the 
physical,  biophysical,  and  biological  factors  that  might  help  explain 
the  mechanisms  of  mammalian  response  to  rapid  variations  in 
environmental  pressure.  Generally,  previous  and  recent  speculations 
have  centered  about  (1)  the  gross  implosion  process,  the  violent  com¬ 
pression  of  the  body  by  an  air  blast  wave,  (2)  the  consequent  changes  in 
volume  of  the  air -containing  organs  due  to  inward  movements  of  the 
thoracic  and  abdominal  walls  and  the  upward  movements  of  the  right 
and  left  diaphragms,  (3)  the  associated  variations  in  internal  pressure 
which  could  be  expected  to  follow  application  of  a  specified  external  pressure 
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load  and,  (4)  the  biological  implications  of  these  several  pressure 
relationships.  4,  8,  12,  22,  27-39 

Tha  present  study  attempts  to  synthesize  these  mechanisms 


in  a  mathematical  model  designed  to  simulate  the  fluid-mechanical  s 

response  of  the  thoraco-abdominal  system  to  sudden  changes  in  f 

environmental  pressure.  Although  the  model  is  simple  in  concept,  \ 

it  was  conceived  to  represent  a  complex  mammalian  organism,  | 

making  it  difficult  to  determine  accurate  parameters  relating  the  •  f 


animal  to  the  model.  However,  even  without  numerical  evaluation, 

;he  model  parameters  can  be  subjected  to  dimensional  and  similarity 

analyses  to  yield  importai.i  relations  between  animal  mass  and  blast  .  ' 

wave  parameters  for  a  given  biological  response.  | 

i 

The  purpoce  of  this  report  is  to  relate  some  of  the  ' 

accomplishments  to  date.  First,  the  mathematical  model  {or  the  \ 

lung  model  for  brevity)  will  be  described.  Second,  a  theoretical  ! 

approach  to  scaling  both  the  body  parameters  of  mammals  and  those  < 

of  the  blast  wave  will  be  presented.  Third,  the  measures  - ~  direct, 
devious  and  arbitrary  --  utilized  to  evaluate  the  parameters  of  the  ? 

lung  model  will  be  described  and  comparisons  will  be  made  between 
pressures  measured  in  the  thorax  of  rabbits  during  blast  and  those 
computed  using  the  model.  Fourth,  examples  of  the  application  of 
the  model  in  interpreting  and  predicting  selected  biological  responses  | 

will  be  presented.  Fifth,  certain  of  the  significant  implications  of  f 

the  work  will  be  pointed  out  and  discussed. 


2.  0 


DESCRIPTION  OF  THE  MODEL 


2.  1  Simplified  Functional  Concept 

An  abnormal  pressure  gradient  is  established  between  gas 
in  the  lungs  and  outside  air  when  the  local  ambient  pressure  changes 
so  fast  that  flow  in  the  airways  cannot  achieve  equalization.  This  is 
particularily  true  when  a  mammal  is  suddenly  engulfed  in  a  blast 
wave  which  is  characterized  by  an  almost  instantaneous  rise  in 
pressure.  T.n  this  case  a  faster  method  of  pressure  equalization  takes 
place  simultaneously  with  flow  in  the  airways;  viz.  ,  reduction  of  lung 
volume  by  the  force  of  the  blast  pressures  acting  on  the  outside  surface 
of  the  body  to  implode  the  thoracic  and  abdominal  walls.  Because  of 
the  inertia  of  the  moving  parts  of  the  body,  overcompensation  will 
occur  and  the  pressures  in  the  lung  will,  at  the  completion  of  the 
implosive  process,  momentarily  oe  high .x  than  those  outside.  Thus, 
an  oscillating  system  is  established  which  is  quickly  damped  because 
of  frictional  losses  of  energy  due  to  tissue  deformation  and  air  motion 
to  and  from,  and  within,  the  lung. 

The  processes  described  above  can  be  stated  mathematically 
and  can  be  evaluated  by  numerical  methods  if  certain  simplifying 
assumptions  are  made.  The  schematic  drawing  in  Fig.  1  illustrates 
the  functional  concepts  of  the  model  which  are  discussed  in  the 
succeeding  paragraphs. 

(1)  In  the  model  a  simple  cavity  of  volume  V  was  used  to 
represent  the  gaseous  volume  subject  to  pressure  changes  upon 
movement  of  the  thoraco-abdominal  structures,  and,  for  the  present, 
let  this  volume  be  that  of  the  lungs  and  airways  which  is  assumed  to 
be  large  compared  to  that  of  gas  in  the  abdomen  in  such  locations 
that  it  could  influence  the  pressure  response  in  the  lungs. 

(2)  Transient  variations  in  volume  of  the  model  cavity  des¬ 
cribed  above  were  ach’eved  through  the  movement  of  two  linear  pistons, 
each  characterized  by  a  mass  M,  an  area  A,  elastic  constant  K,  and 
damping  factor  J.  One  of  the  pistons  represents  the  chest  wall  and 

the  other  that  part  of  the  abdomen  which  moves  with  the  diaphragm.  It 
is  obvious  that  "effective"  values  for  the  above  parameters  should  be 
used.  In  the  case  of  the  mass  and  area  of  the  chest  wall  only  those 
portions  which  move  against  the  lungs  are  involved.  Even  this  is  an 
imprecise  concept  since  the  chest  wall  is  not  flat  ~s  is  the  piston  and 
since  some  parts  of  the  wall,  in  response  to  high  transient  external 
pressures,  will  move  a  ^reat  deal  farther  than  other  parts. 

The  situation  for  the  abdomen  is  even  more  complex.  The 
entire  abdomen  can  be  considered,  for  the  present,  to  be  a  deformable 
"solid"  or  a  viscious  liquid  although  some  gas  is  known  to  be  present. 
When  this  "bag  of  liquid"  is  exposed  to  pressure  or,  its  outer  surface, 
the  "bag"  is  deformed  and  part  of  its  contents  move  with  the  diaphragm 
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Ah 


Model  Equations 


d  x.  dx, 

M,  - -  +  J, - 1 

+  K ,  x  ,  =  A ,  (  P'-P) 

(5) 

1  dt2  '  dt 

d2x,  dx, 

M2 - |  +  J2 - 1 

2  dt  2  2  dt 

+  K2x2=A2(P'-P) 

(6) 

dP  P  dV 

+ 

~!T  "v  "dt" 

1  334xl07  Ak 

h  y 

iP'-PI  cgsu 

(4) 

V 

V  5  V0”X|  A,  - x2  A 2 

(1) 

Fig.  1  Mathematical  model  of  the  thoraco-abdominal 
system  to  simulate  fluid-mechanical  responses  to  rapid 
changes  in  environmental  pressure. 
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into  the  space  originally  occupied  by  the  lungs.  As  with  an  open 
tube  of  toothpaste  under  a  uniform  but  transient  pressure  load  over 
the  cylindrical  surface,  the  greatest  deformation  probably  occurs 
near  the  outlet  (or  diaphragm). 

(3)  In  the  model,  a  simple  orifice  of  equivalent  area  was 
employed  to  represent  the  more  complicated  airways  of  the  mammal. 
Though  it  is  likely  that  the  epiglottis  and  the  vocal  cords  in  reality  act 
to  impede  air  flow  into  and  out  of  the  chest  under  pressure  "loading" 
and  "unloading,  "  the  area  of  the  orifice  was  considered  to  be  such 
that  the  resulting  gas  flow  would  approximate  that  which  actually 
occurs  in  the  animal. 

(4)  The  pressure  and  temperature  of  the  gas  in  the  lungs 
were  assumed  to  be  uniform  over  the  entire  volume  at  any  given  instant 
of  time.  This  assumption  is  valid  if  the  velocity  of  the  boundary  of 

the  lung  is  low  compared  to  the  propagation  velocity  of  pressure  waves 
within  the  lung. 

(5)  It  was  also  assumed  that  the  pressure-volume  changes 

of  the  pulmonary  gas  are  poly  tropic  and,  for  lack  of  definite  information, 
the  polytropic  exponent  was  assumed  to  be  1.  2,  the  average  of  the  ex¬ 
ponents  for  icothermal  and  adiabatic  processes  for  air. 

2.  2  Mathematical  Statement 

If  x  represents  the  inward  displacement  of  a  piston  of  area 
A,  (See  Figure  1)  then  the  cavity  volume  V  is  given  by 

V  =  V  -  x  A  -x  A  (1) 

o  as 

where  V0  is  the  initial  volume  and  the  subscript  1  and  2  ”efer  to  the 
abdominal  and  chest-wall  pistons,  respectively. 

An  expression  relating  the  change  in  cavity  volume  to  a 
change  ir  gaseous  pressure  in  the  cavity,  assuming  a  polycvopic 
process,  is 

dPv/dt  =  -yP/V  dV/dt  (2) 

where  dPv/dt  is  the  time  rate  of  change  of  pressure  resulting  from 
the  rate  of  volume  change  dV/dt,  P  is  gaseous  pressure  in  the  cavity, 
and  y  is  the  polytropic  exponent,  assumed  to  be  1.2. 

Gas  in  the  cavity  also  experiences  a  pressure  change  due  to 
flow  through  the  orifice  of  area  A^.  Since  a  complete  expression  of 
this  phenomenon  is  quite  involved,  an  approximate  relation  obtained 
from  R.  O.  Clark^O  was  u8ed.  The  Clark  equation,  empirically  de¬ 
rived  from  shock  tube  experiments  in  which  cavities  were  filled  by 
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blast-wave  flow  through  orifices,  can  be  written* 

dPf/dt  =  ±  (1.  334  x  107)  (Ah/V)  /|P^rp-p  cg3U  (3) 

The  quantity  dP^dt  is  the  time  rate  of  change  of  pressure  in  the 
cavity  due  to  air  flow,  P'  is  the  outside  pressure,  and  P  is  the 
pressure  in  the  cavity.  The  algebraic  sign  of  the  right  hand  side 
of  the  equation  is  positive  if  the  external  pressure  P'  is  greater 
than  the  internal  pressure  P,  and  negative  otherwise.  Pressures 
are  expressed  in  dynes/cms,  time  in  seconds,  the  area  A.^  in  cm3, 
and  the  volume  V  in  cm3. 

The  net  change  in  the  pressure  of  the  gas  in  the  cavity 
was  assumed  to  be  the  sum  of  Eqs.  (2)  and  (3). 

dP/dt  =  dPv/dt  +  dPf/dt 

=  -yP/V  dV/dt  ±  (1. 334  x  107 )  (Ah/ V)  /]>'  -  PI’  cgsu  (4) 


The  analysis  presented  above  described  cavity  volume  as 
a  function  of  piston  position,  Eq.  (1),  and  gaseous  pressures  in  the 
cavity  as  a  function  of  volume  and  air  flow,  Eq.  (4).  It  now  remains 
to  investigate  piston  motion  as  a  function  of  the  pressure  difference, 
P'  -  P. 

The  piston  experiences  a  driving  force  because  of  the 
difference  in  internal  and  external  pressure  (P*  -  P)  acting  on  the 
piston  area  A.  Opposing  the  driving  force  are  the  forces  due  to 
inertia  M  dax/dta ,  friction  J  dx/dt,  and  elasticity  Kx,  where  M  is 
the  piston  mass,  x  is  the  piston  displacement,  J  is  the  frictional 
force  per  unit  velocity,  and  K  is  the  elastic  force  per  unit  displace¬ 
ment.  Stated  in  equation  form,  the  following  is  obtained: 


M  dax  /dt 8  +  J 

li  l 

dx  /dt  +  K  x  =  A  (P' 

i  ill 

-  P) 

(5) 

M  d3  x  /  dt3  +  J 

S3  3 

dx  /dt  +  K  x  is  A  (P' 
s  a  a  a 

-  P) 

(6) 

The  subscripts  1  and  2  refer  to  the  pistons  representing  the  abdominal 
system  and  the  chest  wall,  respectively.  The  quantities  M,  J,  K, 
and  A,  are  assumed  to  be  constant  for  a  given  mammal.  The  external 
pressure  P'  is  expressed  as  a  function  of  time,  and  the  internal 
pressure  P  is  determined  by  air  flow  and  the  instantaneous  lung  volume. 


*  A  more  accurate  method  for  computing  the  pressure  change  due  to 
air  flow  is  being  developed  and  will  be  described  in  a  subsequent  report. 
The  method  reported  here,  however,  is  satisfactory  for  the  usual  blast 
situation. 
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2.  3  Numerical  Method  of  Solution 


For  convenience  of  solution,  Eqs.  (5)  and  (6)  were  written 
in  terms  of  velocity  v  instead  of  displacement  x  by  the  following 
transformations: 

.  dax/dt3  =  dv/dt 

dx/dt  =  v 

x  =  Cl  v  dt 

oJ 

A  digital  computer  was  used  to  solve  the  model  Eqs.  (1),  (4),  (5), 
and  (6)  simultaneously  for  small  time  steps*  employing  the  following 
approximations: 

(1)  First  derivatives  were  replaced  by  their  finite -difference 
ratios. 

(2)  Functions  were  approximated  by  their  average  values 
over  the  time  interval  assuming  that  the  responses 
were  linear. 

(3)  Integrals  were  approximated  by  their  cumulative  sums, 
summed  over  the  previous  integration  steps. 


*  The  time  steps  used  varied  from  1. 0  to  12  p  sec  depending  on  the 
response  time  of  the  animal  simulated. 
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DIMENSIONAL  ANALYSIS  AND  SCALING 


3.  0 


3 .  1  Model  Equations  in  Dimensionless  Form 

The  model  equations  developed  in  Section  2.2.  are 
summarized  below. 


d^./dt3 

+  dxx  / dt  + 

(K,  /  M.  )  xj  -  (At  /  M,  ) 

{P* 

-  P) 

(7) 

d3xo  ,'dt3 

+  (Ja/M3)  dxg  / dt  + 

(K a/Mg)  Xg  =  (A3/M3) 

(P1 

-  P) 

(8) 

dP/dt  = 

-(yP/V)  riV/dt  ±  (1. 

334  x  iu7)  (Ah/V)  /IT"- 

r-T 

r  cgsu 

(4) 

v  =  vo 

-  xiA,  -  Xg  A3 

(1) 

introducing  the  ambient  air  pressure,  P  ,  and  the  total  body  mass 
u.  3  animal,  m,  various  systems  of  dimensionless  quantities  were 
developed  for  these  model  equations.  The  moat  meaningful  one  is 
presented  below: 


d3  X  i  /  dT3  +  Jj  dX  .  /  dT  4  h^Xi  = 

Ax  (P* 

-  P) 

(9) 

d3  Xa  / dT3  +  Ja  dXa/dT  +  Kg  X3  = 

a3  (P1 

'  P> 

(10) 

dP/dT  =  -(yP/V)  dV/dT  ±  (Ah/V) 

TT 

(ID 

V  =  1  -AiXi  '  A3X3 

(12) 

The  definitions  of  the  underlined  dimensionless  quantities 

(numerics)  are: 

Cavity  Pressure 

P  = 

p/po 

(13) 

External  Pressure 

P>  = 

p'/po 

(14) 

Cavity  Volume 

V  = 

v/v„ 

t  /p  V  u  3  /m 

(15) 

Time 

T  - 

(16) 

Displacement  of  abdominal 

o  o 

piston 

X,  - 

Xl/v^3 

(17) 

Displacement  of  cheat-wail 

piston 

xa  = 

x*/  Vo173 

(18) 

Damping  factor  for  abdominal 
piston 

1^ 

H 

It 

(Ji/Mj)  /m/(F  V17  a) 

(19) 

Damping  factor  for  chest-wall 
pi  ston 

to 

It 

tfa/Ma)  /m/(PoV^  3) 

(20) 

Spring  constant  for  abdominal 
piston 

Kj  = 

K1m/(M1PoVol/  3) 

(21) 
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Spring  constant  for  chest-wall 
piston 

Area  of  abdominal  piston 
Area  of  chest-wall  piston 
Area  of  orifice 

The  dimensional  quantities  are 

P, 

P\ 

P 

o 

Ai  i 

A3  ! 

X1 ' 

X3  . 

V, 

V 

t, 
m, 

Mlt 

Ma  . 

Ah' 

Ji. 

Ja> 

Ki> 

K3i 


Ka  = 

Ks  ni  /  (  hig  P  ■>) 

(22) 

Ai  = 

A1m/(Voa/  3  Mx  ) 

(23) 

Aa  - 

Aam/(VQa/  3  Ma  ) 

(24) 

Ah  = 

1.  334  x  1C'  A  m1/a/(P  V7/  6) 
n  o  o 

cgsu 

(25) 

pressure  in  the  cavity,  or  internal  pressure, 
external  pressure,  a  function  of  time, 
ambient  pressure, 

area  of  abdominal  piston, 

area  of  chest-wall  piston, 

displacement  of  abdominal  piston, 

displacement  of  chest-wall  piston, 

volume  of  the  cavity, 
initial  volume  of  the  cavity, 

time, 

total  body  mass  of  the  animal, 
mass  of  the  abdominal  piston, 

mass  of  the  chest-wall  piston, 

area  of  the  orifice  representing  the  respiratory 
aii  ways  of  the  animal, 

damping  factor  of  the  abdominal  piston,  force  per 
unit  velocity, 

damping  factor  of  the  chest-wal1.  piston,  force  per 
unit  velocity, 

spring  constant  of  tue  abdominal  piston,  force  per 
unit  displacement 

spring  constant  of  the  chest-wall  piston,  force  per 

unit  displacement 


It  should  be  noted  that  the  number  1.  334  "  107  in  Eq.  (25)  has  the 
units  J dynes /s oc  and  hence  caution  must  be  used  with  regard  to  the 
units  of  the  other  quantities  in  the  equation  to  insure  that  a  computed 
A ^  will  be  truly  dimensionless. 

It  will  be  necessary  to  develop  concepts  of  similarity 
among  animals  before  the  significance  of  the  dimensional  analysis 
presented  above  can  be  fully  appreciated. 

3.  2  Animal  Similarity 

It  was  assumed  for  ail  animals  of  a  particular  group  or 
species  that  the  shape  is  the  same  and  that  the  spatial  distributions 
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of  tissue  density,  stretch  (Young's)  modulus,  shear  modulus,  bulk 
modulus,  and  Poisson's  ratio  are  equivalent  (equal  at  corresponding 
points).  Because  equivalence  implies  similarity,  a  single  character¬ 
istic  length  L-,  tissue  density  P0,  stretch  modulus  Y  ,  shear  modulus 
S0,  bulk  modulus  B0,  and  Poisson's  ratio  (tq  are  sufficient  to  specify 
the  animal  for  purposes  of  dimensional  analysis. 


It  then  follows  that 


Ax/Vo3/  3 

- 

Cl 

(26) 

a3/vo3/  3 

r 

c3 

(27) 

Mi  /m 

- 

C3 

(28) 

M3  /  m 

= 

C* 

(29) 

Vc/m 

= 

c5 

(30) 

A.  IVS/  3 
h  o 

= 

Cs 

10  1) 

L/mL'  3 

= 

V“'7 

(32) 

Ki/(Y0D) 

= 

h(S0/\0.  liQ/\0,  <tQ) 

(33) 

K3/iY0L) 

= 

f2  (S0/Y0,  B0/Y0,  <*>) 

(34) 

where  the  f's  are  unspecified  functions  and  the  C's  are  constants. 

(Note  that  the  spring  constants,  Ki  and  Ka.are  inertia-free 
quantities  and  would  therefore  be  independent  of  tissue  density.  ) 
Remembering  that  Y0,  SQ,  BQ,  and  <*Q  are  assumed  to  be  constant 
for  all  animals  in  the  group,  Eqs.  (32),  (33),  and  (34)  can  be 
combined  to  give 

K-./m173  =  C8  (35) 

Ka/ml/3  -  Ce  (  <6) 

The  remaining  animal  paramotcrc  to  be  considered  arc 
the  damping  factors  for  the  abdominal  and  chest-wall  pistons. 

Damping  is  introduced  in  the  model  equations  to  account  for  various 
irreversible  processes  due  to  tissue  deformation  and  air  motion. 

The  experimental  information  on  damping  available  in 
the  physiological  literature  relevant  to  the  pressure  responses  of 
the  chest  and  airways  is  concerned  with  processes  whicti  take 
place  slowly  compared  to  those  resulting  from  exposure  to  air 
blast.  Such  processes  are  slow  enough  to  allow  flow  of  air  to  and 
from  the  lungs,  in  fact,  this  is  an  important  mechanism  involved 
in  the  cited  investigations.  On  the  other  nand,  the  blast-induced 
processes  appear  to  be  so  fast  that  no  significant  flow  can  take  place  in  ‘he 


*  See,  for  example,  Refs.  41,  42,  and  43. 
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airways  although  some  movement  oi  gas  would  be  expected  in  the 
tortuous  passageways  within  the  lung  itself.  Such  flow  would 
be  more  heavily  damped  than  that  in  the  entire  respiratory  tract. 

Thus i  it  is  evident  that  the  damping  processes  involved  in  the  blast 
situation  are  complex  and  cannot  be  rol  ‘ed  directly  to  those  in¬ 
volved  in  conventional  physiological  studies . 

Without  e  precise  under  standing  of  the  damping  processes, 
it  seems  reasonable  to  assume  that  the  "effective"  damping,  or 
damping  ratio,  of  all  species  of  animals  is  the  same,;  i.  e.  ,  for 
all  animals,  the  damping  factor  divided  by  the  critical  damping 
factor  (which  would  just  prohibit  oscillations)  is  the  same.  Ex¬ 
pressed  algebraically, 

i!  =  J/Jc  =  J/  (2/K'M)  =  constant  (37) 

where  o  is  the  damping  ratio,  J  is  the  damping  factor,  Jc  is  the 
critical  damping  factor,  K'  is  the  effective  spring  constant,  and 
M  is  the  mass  being  moved. 

The  coefficient  K1  is  the  sum  of  the  effects  of  (1)  the 
tissue  spring  assumed  to  be  linear  and  described  in  the  model 
equations  as  Ki  and  Ka  and  (2)  the  non-linear  air  spring  resulting 
from  the  fact  that  the  animal's  response  to  blast  is  so  rapid  that  the 
air  initially  in  the  lungs  is  essentially  trapped.  If  it  is  assumed 
that  air  flow  to  and  from  the  lungs  is  not  significant,'  the  pressure- 
volume  relation  can  be  written  as 

PVV  =  PQ  Vo  (38) 

where  P  is  the  pressure  in  the  lung  when  its  volume  is  V.  y  is  the 
polytropic.  exponent,  and  PQ  is  the  pressure  in  the  lung  when  its 
volume  is  V0.  Differentiating  Eq.  (38), 

dP/dV  =  -yP/V  (39) 

The  pressure  differential  uF  is  by  definition  equal  to  -dE/A  and 
the  volume  differential  is  defined  by  the  model  as  -A  dx,  where 
dF  is  the  differential  force  on  the  piston  of  area  A  and  dx  ie  the 
incremental  piston  displacement.  Thus  Eq.  (39)  can  be  written 
in  terms  of  the  non-linear  spring  "constant"  as 

dV/dx  =  -yAsP/V  (40) 

The  variable  volume  V  can  be  eTminated  from  Eq.  (40)  by  in¬ 
troducing  Eq.  (38). 

Vy 

cLF/dx  =  -(yA2P/V0)  (P/PQ)  (41) 

*  Experimental  evidence  pertaining  to  air  flow  during  the  blast 
experience  can  be  found  in  Refe.  4  and  44. 
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The  effective  spring  coefficient  can  now  be  expressed  as 


K1  '  K  -  dF/dx 

=  K  +  ( \Ar  P /  V  )  (P/P  )1/Y  (42) 

o  o 

so  that  lhqs.  (42)  «ind  (37)  can  be  combined  to  give 

J  =  2cM17  2  [K  4  (vA3P/ V  )  (P/P  )l/v  YJ  2  (4i) 

o  o 

Matting  use  of  the  scaling  relations  in  Eqs.  (26)  -  (36),  Eq.  (43)  can 
be  written  for  tne  two  pistons 

J  /  m?/  3  =  2«  [C  C  +  yC2  C  CP7  3  P  (P/P  )’'/Vj17  2  (44) 

l  3  8  13°  O 

J  /  mr/  3  =  2a  [  C  C  +  v  C2  C  C17  3P  (P/P  2  (46) 

In  the  absence  of  direct  evidence  it  was  assumed  that  the  damping 
ratio  has  a  constant  value  a0  for  all  animal  species  when  both  P  and 
P  have  a  value  of,  say,  1.013  x  10  dynes/  cm*  (14.  7  psi).  Eqs.  (44) 
and  (45)  can  then  be  written 

J  /  m!/  3  =  2o  f  C  C  +  yC2C  C17  3  (1 . 013  x  10e)  Ji7  3  cgsu  (4o) 

1  o  3  a  1  i  3  s  ' 

J  /ma'  3  =  2 a  [  C  C  +  y  C2C  C1'  3(1. 013  x  106)]17  2  cgsu  (47) 

2  O  4  9  2  4  15 

or 

J  /  m8'  3  =■  C 

1  10  (48) 

J  / m2/  3  =  C  (49) 

r  1 1 


It  is  of  interest  to  note  for  a  given  animal  (  J  and  m  both  constant) 
the  change  in  the  damping  ratio  a  caused  by  changes  in  P  and  P  ,  an 
increase  in  P  results  in  an  increase  in  the  damping  ratio  but  an  in¬ 
crease  in  P  lowers  it.  See  Eq.  (43).  During  oscillation  the  average 
internal  pressure  is  about  the  same  as  the  external;  thus,  for  exposures 
at  a  given  ambient  pressure,  the  thoraco-abdominal  system  would 
vibrate  more  freely  in  response  to  a  blast  wive  of  "long"  duration 
than  it  would  to  One  of  "short"  duration  --  the  durations  being  "long" 
or  "short"  relative  to  the  period  ot  oscillation. 

Similarly,  "long"  -  duration  blast  waves  of  small  magnitude  would 
prodi.ee  osr  llations  with  greater  damping  than  would  corresponding 
waves  of  larger  magnitude. 

3.  3  Environmental-Pressure  Scaling  for  Similar  Animals  of 
Different  Masses 

The  model -equation  numerics  defined  in  Eqs.  (13)  -  (25)  can 
be  sir  plified  by  using  the  similarity  relations  in  Eqs.  (26)  -  (36),  (48), 

and  (<i  ))  to  eliminate  M  ,  M  ,  K  ,  K  ,  J  ,  J  ,  A  ,  A  ,  A,  ,  and  V 

1  21  I3h  o 


12 


as  scaling  parameters.  For  completeness,  Eqs.  (13)  and  (14) 
are  repeated. 


P  =  (P/P  )  (13) 

—  o 

P'  =  (P'/P0)  (14) 

V  =  (V/m)  (Cjl)  (44) 

T  =  (tP1''  3/mx/  3)  {Cl/  G)  (45) 

—  o  b 

X  =  (x  /m”  3)  (C"1/  3)  (46) 

—  1  1  O 

X  =  (x  /in”  3)  (C'U  3)  (47) 

—  3  2  B 

J  =  (l/P17  “)  (C_1  C_1/  3  C  )  (48) 

—  1  O  3  E  10 

J  =  ( 1  / Pq  s )  (C-1  C’U9C)  (49) 

K  =  (1/P  )  (C"1  C_1/  3  C  )  (50) 

*—  1  O  3  6  Q 

K  =  (1/P J  (C~l  C'U  3  C  )  (51) 

A  =  (C  C'1)  (52) 

“I  13 

A  =  (C  C'1)  (53) 

Z  3  3 

A.  =  (1/P)  (  1 .  334  x  1 07  C  C'1'  3)  cgsu  (54) 

n  O  6  6 


Note  that  the  animal-parameter  nun«cricr>  defined  by  Eqs.  (52) 
and  (53)  have  constant  values  (independent  oi  animal  mass  and  ambient 
pressure),  the  remaining  animal  numerics,  Eqs.  (48)  -  (51)  and  (54), 
are  independent  of  animal  mass  but  do  depend  on  ambient  pressure. 

The  scaling  relations  involving  changes  in  animal  mass,  which  are 
derived  below,  are  based  on  the  condition  that  the  animal-parameter 
numerics  are  independent  of  animal  size  anrl  thus  on  the  validity  of  the 
similarity  relations  derived  in  the  last  section.  However,  to  derive 
the  scaling  relations  involving  changes  in  the  ambient  pres  sure  it  is 
necessary  to  neglect  the  effect  of  these  changes  on  the  spring-constant, 
damping,  and  the  oi  ifice -are.-  numerics  mentioned  in  tile  previous 
paragraph.  Having  made  these  approximations,  all  of  the  coefficients 
of  the  terms  m  Eqs.  (9'  -  ( IP)  become  constant  (independent  of  ambient 
pressure  and  animal  r  iss)  so  that  if  we  consider  two  cases  (differing 
in  animal  mass  or  am'  if  t  pressure,  or  both)  where  the  initial  scaled 
piston  positions  (X  1  ana  and  piston  velocities  (dXi  / dT  and  dX2  /dT) 
are  equal?'  the  corresponding  functions  P  (T),  V(T  ),  Xi(T),  and3fa/(T) 
will  all  be  equal  between  the  two  cases  provided  only  that  the  input 
functions,  P  1  (  T )  ,  are  equal,  i.  e.  ,  if  the  scaled  initial  conditions  are  equal  and 

*  Note  that  these  conditions  also  require  the  initial  scaled  volume  (V) 
and  rate  of  change  of  volume  (dV  /dT)  to  be  equal.  “ 
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then 


where  the  subscripts  1  and  2  relei  to  the  first  and  second  case, 
respectively. 


According  to  Eq.  (5o),  the  internal-pressure  numerics 
as  functions  of  the  time  numerics  are  the  same  for  two  exposure 
siutations,  where  the  animal  masses  and  /or  the  ambient  pressures 
are  different,  provided  the  scaled  pressure-time  input  parameters 
are  equal  as  stated  in  Eq.  (55).  Examples  of  the  application  of 
these  equations  to  problems  in  blast  biology  are  presented  in 
Sec.  5.0. 


(55) 


(5b) 
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•5.  0 


EVALUATION  OF  THE  ANIMAL  PARAMETERS  USED  IN 


THE  MODEL  EQUATIONS 


4.  1  General  Remarks 

An  attempt  was  made  to  establish  reasonable  estimates 
of  the  model  parameters  by  various  methods  described  below. 

The  set  of  estimated  parameters  were  then  tested  by  using  them 
in  the  model  equations  to  compute  the  thoracic  pressure  patterns 
for  rabbits  exposed  to  shock-tube  blast  waves  and  comparing  the 
results  with  those  obtained  experimentally.  Adjustments  were 
then  made  in  one  or  more  of  the  parameters  and  the  comparison 
was  repeated.  Several  such  trials  were  made  until  the  computed 
pressure-time  patterns  corresponded  reasonably  well  with  those 
measured.  A  summary  of  the  results  is  presented  in  Sec.  4.  7. 

4.  2  Lung  Volumes  and  Densities 

45 

Results  of  a  study  by  Crosfill  and  Widdicombe  in>  luded 
measurements  of  functional  residual  capacity,  tidal  volume,  and 
lung  mass  for  seven  species  of  animals.  The  report  also  contains 
data  for  man  taken  from  other  source?.  46,  47 

According  to  the  authors  the  techniques  used  to  measure 
the  functional  residual  capacities  of  the  lungs  tended  to  under¬ 
estimate  the  true  volumes.  However,  since  the  same  methods  were 
used  for  all  species,  the  results  are  valuable  for  purposes  of 
comparison.  Average  gaseous  volume  of  the  lungs  per  unit  body 
mass,  and  average  lung  density  are  plotted  in  Fig.  2  as  functions 
of  body  mass.  (These  quantities  are  defined  in  the  figure.  )  Note 
that  the  data  for  the  mouse,  rat,  guinea  pig,  and  rabbit  are  sig¬ 
nificantly  different  from  those  for  the  monkey,  cat,  dog,  and  man. 
The  average  of  the  average  lung  volumes  per  unit  body  mass  is 
9.  08  cm3/kg  for  the  smaller  species  and  29.  8  cm3 /kg  for  the 
larger.  However,  the  average  lung  density  for  the  small  species 
is  almost  twice  as  great  as  for  the  large,  0.  367  and  0.  194  g/cm3, 
respectively. 

While  the  data  in  Fig.  2  are  useful  in  estimating  average 
lung  volumes  for  animals  of  either  the  small  or  the  large  species, 
it  is  evident  that  scaling  between  the  two  groups  by  body  mass 
would  result  in  gross  errors.  Average  lung  densities  are  used  in  the 
next  section  to  help  estimate  piston  areas  and  masses. 

4.  3  Geometric  Model  of  the  Thorax 

in  the  last  section  estimates  were  made  for  the  average 
gaseous  volume  of  the  lungs  and  average  lung  density,  making  it 
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V/m,  cm3 /kg 


E 

o> 

i«C 


Body  Mass(m),  kg 


V/m*  (FRC+*i/2TV)/m 
V  >  Avg  gaseous  volume  of  lungs 
m  >  Body  mass 
FRC  s  Functional  residual  capacity 
TV  1  Tidal  volume 
M  >  Mass  of  lungs 


p=  M/(Vtv) 
p  *  Avg  density  of  lungs 
v  *  Volume  of  lung  tissue  * 

*  D  snsity  of  tuna  tissue  assumed  to 
be  I.Og/cm3 


Fig.  2  Average  lung  volume  per  body  mass  and  average  lung 
density  as  functions  of  body  mass  for  seven  species  of 
animals  and  man.  The  plots  were  prepared  using  data 
in  Ref  45.  Lung  mass  as  a  percent  of  body  mass, 
lOOM/m  =  IOOtt  (  V/m)  /  (1  -  ?)  ,  is  0.  526%  for  the 
small  animals  and  0,  717%  for  the  large  ones. 
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possible  to  compute  the  average  total  lung  volume  (gas  plus  tissue). 

In  this  section  an  attempt  will  be  made  to  relate  piston  areas  and 
masses  with  lung  volume  --  the  relation  being  made  more  realistic 
by  idealizing  the  thorax  as  a  regular  geometrical  shape. 

48 

In  a  meticulous  work,  Eycleshymer  and  Shoemaker 
illustrated  and  described  cross -s ections  of  human  cadavers  at  in* 
tervals  of  2  to  3  cm.  The  sections  passing  through  the  thorax  were 
used  hero  to  determine  the  cross-sectional  areas  of  the  lungs,  the 
circuinl  rential  lengths  of  contact  of  the  lungs  with  the  rib  cage,  and 
the  cross-sectional  areas  of  the  rib  cage  in  contact  with  the  lungs. 

Note  that  the  last  two  measurement,  involve  only  those  portions  of 
the  chest  wall  which  could  be  expected  to  move  enough  under  dynamic 
loading  to  reduce  the  volume  of  the  lungs  significantly. 

If  the  cross-sectional  areas  of  the  lungs  are  assumed  to  be 
circular,  equivalent  diameters  can  be  computed  for  the  measured 
areas  and  also  for  the  measured  lengths  of  contact  with  the  rib  cage. 

Fig.  3  is  a  plot  of  these  equivalent  diameters  for  the  9  cross-sections 
passing  through  the  lungs.  Note  that  the  equivalent  diameters  based 
on  each  of  the  criteria  increase  with  height  above  the  lower  lobe  .rom 
zero  to  oroximately  the  same  maximum  values  in  the  region  of 
the  diaphiagm.  The  maximum  area  of  the  lung  cross-section  in  this 
region  (230  cm3)  can  be  taken  here  as  an  approximation  of  the  effect¬ 
ive  area  of  the  piston  representing  the  abdomen.  The  total  volume 
of  the  lungs,  gas  plus  tissue,  was  determined  by  summing  the  products 
of  cross-sectional  areas  and  the  distances  separating  the  cross  sections. 
Similarly,  the  products  of  circumferential  lengths  and  separation 
distances  were  summed  to  give  tne  total  surface  area  of  contact  of 
lungs  with  rib  cage.  These  proved  to  be  2960  cm3  and  704  cnT" , 
respectively. 

By  trial  and  error  it  was  found  that  a  2960  cr.  3  right-circular 
cylinder  whose  height  is  3/4  of  its  diameter  has  the  desired  area 
relations  withir  the  accuracy  of  the  measurements.  The  following 
tabulation  illustrates  the  point: 

measured  using  the  computed  for 

cross-sections  right-circular 

cylinder 


V+,  cm3 

2960 

2960 

Ai ,  cma 

230 

230 

As  ,  cm2 

704 

691 

V  i6  the  volume  of  the  gas  and  tissue  in  the  lungs:  Ai , 
the  area  of  the  abdominal  piston,  is  the  area  of  the  base  of  the  cylinder; 
ani  Aa  ,  the  area  of  the  chest-wail  piston,  is  the  area  of  the  curved 
surface  of  the  cylinder.  The  diameter  and  the  height  of  the 


Height  Above  the  Lower  Lobe  of  the  Lungs,  cm 


Diameter  of  the  Lung  Cavity,  cm,  Assuming  Circu¬ 
lar  Cross- Sections  Based  on  the  Measured  - 

(1)  •  Area  of  the  lungs 

(2)  a  Length  of  contact  of  the  lungs  with  the 

rib  cage 

Fig.  3  The  design  o£  a  geometric  model  of  the  thorax  using 
cro££-ssc*ioncl  del ici  £ o r  hu’n 3.n  c ? r s 

in  Ref  48. 
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right-circular  cylinder  are  shown  in  Fig.  3  as  dashed  lines. 


The  volume  of  the  rib-cage  tissue  adjoining  the  lungs, 
measured  from  the  cross-sections  and  their  separations,  was 
2744  cm3.  Assuming  that  the  tissue  had  a  density  of  1  gm/cm  , 
the  mass,  M3,  of  the  piston  representing  the  chest  wall  yjjas  estimated 
to  be  2744g.  For  tne  lack  of  more  definitive  information,  the 
abdominal  piston  was  assumed  to  have. the  came  rr.ts:  per  unit 
area  as  the  chast-wall  piston. 


The  relations  derived  above  can  be  expressed  as  scaling 
equations  to  help  estimate  credible  relationships  between  the 
gaseous  and  total  volumes  of  the  lungs  and  the  areas  and  masses 
of  the  two  pistons  taking  part  in  the  implosive  process . 


giving 


> 

II 

1.117  V+a/a 

(57) 

As 

3.  352  V  a'  3 

+ 

(58) 

Mj 

0.  3090  V  p 

T  t 

(59) 

Ma 

0.  9270  V+Pt 

(60) 

V+  = 

V  /  ( 1  -  v/  Pt) 

(61) 

Ai/Mi  = 

Aa/Ma  =  3.  6 1  5/(  3  p  ) 

(62) 

The  symbols  used  are  defined  blow. 


A; 

A3 

Mx 

Ma 


V 

F 

n 

"t 


area  of  abdominal  piston 

area  of  chest-wail  piston 

mass  of  abdominal  piston 

mass  of  chest-wall  piston 

average  total  volume  of  lung  cavity  (gas 

plus  tissue) 

average  gaseous  volume  of  lungs,  see  Fig. 
average  density  cf  lungs,  see  Fig.  2 
tiw sue  dsnsityi  2.^ s ur^isd  tc 


Note  that  the  areas  and  masses  ir.  Eqs.  (57)  -  (62)  are 
scaled  with  the  average  total  lung  volume  (  a  function  of  average 
gaseous  volume,  average  lung  density,  and  tissue  density)  rather 
than  with  body  mass.  This  procedure  facilitates  use  of  the 
scaling  data  in  Fig.  2  for  the  two  distinct  groups  of  species. 


4.  4  Spring  Constants 

41  42 

Velocity  resonant  irequencies  have  been  measured  fer  cats,  ’ 

^  See  discussions  in  Sec3.  2,  1  and  6.0. 
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and  dogs  by  introducing  sinusoidal  waves  into  the  r : — v  a  y  s  via 
the  trachea  or  by  exposing  the  whole  animal  to  oscillating  pressures 
ir.  a  closed  chamber  with  a  tracheal  tube  connected  to  measuring 
devices  outside  the  chamber.  In  either  case  the  forcing  fre¬ 
quencies  were  varied  to  determine  that  frequency  which  would 
produce  either  maximum  chest-wall  velocity  or  maximum  air 
velocity.  Tne  velocity  i  esonant  frequency  is  defined  by  the  follow¬ 
ing  equation: 

fv  =  /K7M/(2n)  (63) 

where  fy  is  the  frequency  producing  the  maximum  velocity  re¬ 
sponse,  K  is  the  spring  constant  with  units  of  force  per  unit 
displacement,  and  M  is  the  mass  moved. 

The  response  due  to  motion  of  the  chest  wall  was  not 
separated  from  that  due  to  motion  of  the  diaphragm  and  abdomen, 
i.  e.  ,  a  single  resonant  frequency  was  determined  for  each  animal. 

Thus,  to  apply  these  data  to  the  model  parameters,  M  in  Eq.  (63) 
would  represent  the  mass  of  either  of  the  pistons  and  K  the  corres¬ 
ponding  spring  constant.  Since  mass  per  unit  area  was  assumed  to 
be  the  same  for  the  two  pistons  (see  Sec.  4.  3),  the  spring  constant 
per  unit  area  would  also  be  the  same  if  the  pistons  have  the  same 
resonant  frequency. 

To  make  the  resonant  frequency  data  for  cats  and  dogs 
more  generally  applicable,  it  is  necessary  to  determine  a  mass¬ 
scaling  relation.  According  to  Eqs.  (35)  and  (36),  the  spring  constant 
varies  as  the  cube  root  of  the  animal  mass,  and  referring  to 
Eqs.  (28)  and  (29)  the  vibrating  masses  vary  directly  as  the  body 
mass.  Applying  these  relations  to  Eq.  (63),  the  following  is 
obtained 

fyml/  3  =  /C 77CT  /  (2*)  ^  /C TclnZn)  =  Cia  (64) 

where  m  is  the  body  mass  and  C  is  a  constant. 

Eq.  (64)  was  used  to  normalize  the  resonant  frequency 
data  from  references  41,  42,  and  43.  The  results  are  shown  below: 


Animal 

f 

v  -l 
sec 

m 

kg 

f  3 

1  m 

V 

1/  3  / 

g  /sec 

Ref, 

cat 

9.  6 

2.  50 

130 

41 

cat 

11.5 

2.  76 

161 

42 

dog 

5.  4 

12. 

124 

43 

Avg.  138 


i 
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rT  Making  use  of  Eqs.  (63)  and  (64),  the  following  scaling 

w  equations  are  obtained  by  equating  fvrrr  3  to  the  average  experimental 

T  value  of  138  3  sec'1  : 

'  r  Ki  /  Mi  =  Ks  /  -  751,800  /  m2/  3  cgsu  (65) 

f  "~T  « 

j 

!  The  subscripts  1  and  2  refer  to  the  abdominal  and  chest-wall 

i  -  '  pistons,  respectively. 

I 

I  4.  5  Effective  Orifice  Area 


Results  of  rapid-decompression  tests  on  a  83-kg  man 
•were  reported  by  Euft  and  Bancroft. 49  Using  an  analytical 
technique  developed  by  Haber  and  Clamaim^O  the  authors 
evaluated  an  average  effective  orifice  diameter  of  0.  478  cm. 

By  applying  Eq.  (30)  and  (31),  the  following  scaling  relation  is 
obtained  for  the  effective  orifice  area. 

A ^  -  9.  43  x  10  mS/  3  cgsu  (66) 

4.  6  Damping  Factors 

Experiments  in  blast  biology ^ *  produced  records  of 
pressures  in  the  thorax  of  animals  exposed  to  blast  waves  generated 
by  a  shock  tube.  Two  such  records  are  shown  in  Figs.  4  and  5. 

In  each  of  the  experiments  a  small  pressure  transducer  was 
inserted  in  the  esophagus  of  a  1.  8-kg  rabbit  down  to  the  level 
of  the  heart.  The  animal  was  then  placed  in  a  cage  bolted  side-on 
against  the  end-plate  closing  the  shock  tube.  In  each  case,  the 
upstream  side  of  the  animal  "felt"  the  incident  shock  of  about 
7-psi  overpressure  followed  by  the  reflected  shock  approximately 
0.  5  msec  later  --  the  travel  time  of  the  shockwave  across  the 
animal  to  the  reflecting  surface  and  back  again.  However,  the 
side  of  the  animal  against  the  enu-plate  "felt"  the  incident  and  the 
reflected  shocks  almost  simultaneously.  The  time  step  used 
in  the  computations  was  the  average  of  those  for  the  limiting  cases 
mentioned  above;  viz,  0.25  msec,  representing  the  time  separation 
of  the  incident  and  reflected  shocks  at  a  distance  from  the  end- 
plate  equal  to  that  of  the  cc-nter  of  the  thorax. 

The  internal- pi  es sure  records  described  above  are 
characterized  by  damped  oscillations  and  were  used  to  estimate 
damping  factors.  For  analytical  reasons  it  was  necessary  to  assume 
a  linear-spring  system  even  though  the  "air  spring"  of  the  lung 
model  makes  the  system  non-linear.* 


*  See  discussion  in  Section  3.  2 
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Overpressure,  psi 


(67) 


According  to  linear  theory  the  formula  for  the  damping  relation  is 

-1/  s 

a  =  In  (x  /x  )  [l+lns(x/x  )  /  (4tt3  )  ] 
n  n+i  n  ni-i 

2tt 

where. g  is  the  ratio  of  the  damping  factor  to  the  critical  damping 
factor,  x  and  *n+,  are  the  amplitudes  of  two  consecutive  vibrat¬ 
ions,  ancPln  indicates  the  logarithm  to  the  base  e. 

Information  available  at  this  time  will  not  allow  the  effects 
of  the  two  pistons  (abdominal  and  chest-wall)  to  be  distinguished 
from  each  other.  Thus  the  simple  assumption  was  made  that  the 
pistons  move  as  one;  i.  e.  ,  their  respective  displacements  are 
identical  at  any  point  in  time  provided,  of  course,  that  their 
initial  displacements  and  velocities  are  the  same.  Examination 
of  the  model  equations,  (7)  and  (8),  shows  that  the  conditions  for 
the  two  pistons  to  move  as  onev  are  that  K1/M1  =  Ks/Ma, 

Ai  /  Mj  =  Aa  /Ma  ,  and  Ji  /M-^  =  Ja  Ma  .  The  first  two  conditions 
were  dealt  with  in  previous  sections  (4.  4  and  4.  3',,  and  the  third 
one  is  recognized  in  the  material  that  follows. 

In  order  to  use  Eq.  (67)  it  is  necessary  to  express  the 
displacement  amplitudes  in  terms  of  peak  pressures  which  were 
measured.  Neglecting  air  flow,  the  polytropic  pressure -volume 
changes  of  the  lung  can  be  expressed  as 

x  =  V  / A  [1  -  (P  /P)  V ] 

where  x  is  the  displacement  of  each  of  the  pistons  with  total  area 
A  when  the  cavity  pressure  (absolute)  is  P,  V0  is  the  initial  volume 
of  the  cavity  when  the  internal  pressure  ii  PD  and  the  piston 
displacement  is  zero,  and  y  is  the  polytropic  exponent  assumed 
to  be  1.2. 


Ddid  for  the  firsi  two  peaks  of  the  measured  internal 
pressure  records  shown  in  Figures  4  and  5  were  used  in  Eqs. 

(67)  and  (68)  to  evaluate  damping  ratios.  The  succeeding  peaks 
could  not  be  used  since  the  blast  pressures  were  decaying 
appreciably  afte:  the  second  peak.  Damping  ratios  of  0.  149  and 
0.  130  were  evaluated  for  the  data  in  Figs.  4  and  5,  respectively, 
while  the  internal  pressure  was  oscillating  about  an  external  over¬ 
pressure  of  approximately  20  psi  iji  each  case. 

The  average  value  (0.  1395)  of  the  damping  ratios  mentioned 
above  were  used  in  Eq.  (43)  along  with  animal  parameters  previously 
estimated;  the  total  damping  factor  thus  determined  for  the  1. 8-kg 
rabbits  was  37,100  dynes -sec/ cm.  As  mentioned  above,  damping 

*  The  single  piston  equivalent  to  the  two -piston  system  has 
parameters  equal  to  the  sum  of  those  for  the  two  pistons;  i.  e.  , 
Mi+Ma,  Ax+  A3 1  etc. 


(68) 
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was  apportioned  to  each  uf  the  pistons  according  to  its  mass, 
resulting  in  the  following  scaling  relation. 

Ji/Mi  =  J3/M3  =  14  1  70/m1/f  3  cgsu  (6?) 

where  Jj  and  Jy  are  damping  factors  for  the  abdominal  and 
chest-wall  pistons  of  mass  "M.  and  M.,,  respectively  and  m  is  the 
animal  body  mass. 

4.  7  Summary  of  Estimated  and  Adjusted  Parameters 

It  was  shown  in  previous  sections  that  the  body-mass 
scaling  of  most  of  the  animal  parameters  depends  on  the  species 
type.  Evidence  was  presented  of  at  least  two  types:  the  mouse, 
rat,  guinea  pig,  and  rabbit  being  members  of  the  small-species 
type  and  the  monkey,  cat,  dog,  and  nan,  members  of  the  large- 
species  type. 

Estimated  model  parameters  scaled  to  an  idealized 
1-kg  animal  are  listed  in  the  first  two  columns  of  Table  1  for 
the  two  types  of  species.  The  parameter  which  differs  greatest 
between  the  two  types  is  that  of  average  lung  volume,  the  large- 
species  type  having  an  average  volume  more  than  three  times  as 
great  as  that  for  the  other.  It  is  not  significant  that,  the  estimated 
orifice  areas  are  the  same  for  the  two  types  since  they  were  both 
scaled  according  to  body  mass  from  data  obtained  for  man,  ' 
additional  information  will  probably  indicate  they  are  actually 
different. 


The  third  column  of  Table  1  lists  the  model  parameters 
scaled  to  a  1  8-kg  small-species  animal.  The  last  column  of  the 
table  contains  the  "adjusted"  values  of  the  parameters  obtained 
by  attempts  to  match  the  output  of  the  model  to  experimental  records 
of  pressure  in  the  thorax  cf  1.  8-kg  rabbits  during  blast  experiences. 
(See  Sec.  4  6)  The  pressure-time  curves  marked  "computed” 
in  Figs.  4  and  5  were  obtained  using  these  adjusted  parameters. 
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Tabic  1  SUMMARY  OF  ESTIMATED  AND  ADJUSTED  MODEE  PARAMETERS 


Estimated  for  1-kg  Animall 

1 . 8 -kg  Rabbit  j 

Reference 

"Small" 

Species 

"Large" 

Species 

Estimated 

Adjusted 

V,  cm3 

9.  08 

29.  8 

16.  3 

10.  3 

Fig  2 

Ai«  cm3 

o.  59 

12.4 

9.  74 

8.  78 

Eqs  (57)  (61) 

Aa  >  cm2 

19.  8 

37.2 

29.  2 

26.  3 

Eqs  (98)  (bl) 

Ml  g 

4.  43 

11.4 

7.  9b 

7.  04 

Eqs  (59)  (61) 

In  g 

13.  3 

34.  3 

23.  9 

21. 1 

Eqs  (bO)  (61) 

Ki  dynes  / cm 

33,  300 

85,  700 

40,400 

39, 900 

Eq  (65) 

K3  dynes /era 

100,  000 

258.  000 

121, 000 

12  0,  000 

Eq  (OS) 

Ah  cm2 

0. 00943 

0.  00943 

0. 0140 

0.  0140 

Eq  (6b) 

JT  ,  dynes-sec/cm 

6,  280 

lb,  200 

9.  275 

6,  570 

Eq  (o9; 

Ja  ,  dynes  -sec/  cm 

18, 800 

48,  600 

27,  825 

19, 700 

Eq  (b9) 

Definitions 

V 

A; 

Aj 

M; 

Ma 

Ki 

Kj 


Ji 

•fa 


Average  gaseous  volume  of  the  lungs 

Area  ul  piston  representing  the  abdomen 

Area  of  piston  representing  the  chest  wall 

Mass  of  piston  representing  the  abdomen 

Mass  of  piston  representing  the  chest  wall 

Spring  constant  for  the  piston  representing  the  abdomen 

Spring  constant  for  the  piston  representing  the 

c lie st  wall 

Effective  area  ct  the  orifice  representing  the  airways 

Damping  factor  for  the  piston  representing  the  abdomen 
Damping  factor  for  the  piston  representing  tne 
chest  wall 
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APPLICATIONS  TC  PROBLEMS  IN  BLAST  BIOLOGY 


5.  0 


5.  1  Restatement  of  the  Scaling  Equations 

The  scaling  equations  derived  in  Sec.  3.3  for  similar 
animals,  (35)  and  (56),  apply  to  both  the  external  and  internal  pressure 
waves  as  functions  of  time.  For  convenience  in  making  specific 
applications,  the  equations  are  written  below  in  terms  of  characteristic 
pressures  and  times  used  to  compare  pressure  waves  of  the  same  shape. 
(Pressure  waves  are  considered  to  be  of  the  same  shape  if  one  can  be 
derived  from  the  other  by  multiplying  the  pressures  and  the  corresponding 
time9  by  constant  factors,  a  pressure  factor  and  a  time  factor.) 


P'  /  (P)  =  P’/  (P  ) 

z  o  a  l  o  i 

t  (P  )i;  3  /  m  1/3  =  t  (P)l/  a  /  my  a 
a  o  j  a  i  o  i  i 

P  /(P)  =  P  /  (P) 

a  o  3  l  o  a 


(69) 

(70) 

(71) 


where  the  subscripts  1  and  2  refer  to  two  blast  situations  where  the 
blast  parameters  and/or  the  animal  body  masses  are  different.  The 
parameter  definitions  are: 

P'  :  Any  characteristic  pressure  of  blast  waves  of 

similar  shape.  Usually  the  peak  pressure  is  used 
to  compare  two  different  blast  waves. 

P  :  Ambient  pressure,  also  assumed  to  be  the  initial 

°  gaseous  pressure  in  the  lungs. 

P  :  Any  characteristic  pressure  of  the  pressure  pulse 

occurring  in  the  lungs.  The  peak  lung  pressure  is 
used  in  this  study. 

t  :  Any  cha rac te r i stir  time  of  blast  waves  of  similar 

shape,  or  ol  the  internal  pressure  waves;  e.  g.  , 
duration  of  the  biast  wave,  time  to  reach  maximum 
pressure  in  the  animal  lungs,  etc. 

m  ;  Body  mass  of  similar  animals.  See  Sec.  3 . 2 

5 .  2  Index  of  Biological  Respon s e 

In  the  design  of  the  model  the  input  numerics  J_i ,  £  a  .  , 

K3,  Ai  ,  Aa  ,  A,  ,  and  P  (T)  were  assumed  to  be  ail  the  significant 
parameters  needed  to  describe  the  gross  response  of  the  thoraco¬ 
abdominal  system  tc.  rapid  changes  in  environmental  pressure.  Having 
made  this  assumption,  it  follows  from  dimensional  analysis  that  any 
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nondimensional  index  of  biological  response  R(  e.  g.  ,  percent  mortality, 
percent  increase  in  lung  mass,  etc.)  i6  a  function  only  of  the  numerics 
listed.  Thus,  letting  the  bracket  subscripts  1  and  Z  refer  to  two  blast 
situations  where  the  blast  parameters  and/or  the  animal  body  masses 
are  different,  it  can  be  stated  that 


if  [J_i  !i  =  [Ji  ]a  *  [Ja  j  1  -  [ J?  Ja  .  I.Ki  ji  -  [Ki  la  , 

[K?  Ji  =  [Kr?  la  ,  ■"A  j  li  =  [A i  J3,  '.Aj  Jr  =  [As  Js  ,  and 
[A  hli  =  [Apia,  as  was  assumed  in  Sec.  3.  3  in  deriving 

scaling  equations  (55)  and  (56)  for  similar  animals;  and  if  LP'(T)  Ji  = 

[p1  ( T)  ]a;  1.  c.  ,  if  Eqs .  (b9)  and  (70)  are  satisfied. 

Note  that  the  above  analysis  does  not  indicate  specific  mechan¬ 
isms  but  only  the  relations  between  the  sigmlicanl  physical  parameters 
influencing  biological  response;  however,  mechanisms  may  be  postulated 
by  interpretation  of  Eqs.  (69)  -  (7Z).  Subject  to  the  restrictions  previously 
discussed,  two  mammals  exposed  to  different  blast  situations  will  have 
equal  biological  response,  Eq.  (72),  if  the  conditions  of  exposure  satisfy 
Eqs.  (69)  and  (70).  Thus,  referring  to  Eq.  (71),  the  ratio  of  the  peak 
internal  pressure  to  the  ambient  pressure  (P/PG)  could  be  the  critical 
factor  to  the  animal;  however,  it  should  be  noted  that  this  ratio  is 
approximately  proportional  to  the  r£,tio  of  the  initial  lung  volume  to 
the  minimun  volume  (V0/V),  a  measure  of  tissue  distortion.  Another 
critical  factor  could  be  the  ratio  AP/P„  where  AP  is  the  difference 
between  the  peak  internal  pressure  and  the  simultaneous  external  pressure, 
and  Pc  is  the  ambient  pressure.  The.  factors  which  are  critical  in  pro¬ 
ducing  biological  damage  from  air  blast  may  be  determined  more 
precisely  in  future  investigations. 

5.  3  Long -duration  Blast  Waves 

Six  species  of  animals  (mouse,  rat,  guinea  pig,  rabbit,  dog, 
arid  goat)  were  used  in  ouk  shuck,  tube  study*0  ana  tv.-c  species  (hamster 
and  cat)  in  another-**  to  determine  biological  tolerance  to  sharp-rising 
blast  waves  of  "long"  duration  (400  msec).  Except  for  the  dogs  and  goats 
which  were  secured  with  harness  and  straps,  the  animals  were  exposed 
in  wlJe-mesh  cages  against  the  reflecting  surface  at  the  tnd  of  the  shock 
tube.  Thus,  the  upstream  side  of  the  creature  felt  the  incident  wave 
first  ana  a  short  time  later  the  reflected  wave.  The  length  of  this  time- 
step  would  be  approximately  proportional  to  cue  body  width  of  the  animal 
or  to  the  cube  root  of  the  body  mass.  Thus,  by  Eq.  (  70),  the  effect  of 
the  time  step  would  be  identical  for  (similar)  animals  of  all  body*  masses 

Examination  of  the  measured  pressure-time  records^  reveals 
that  the  blast  pressures  were  approximately  constant  during  the  time 
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cf  interest,  i.  e.  ,  the  first  few  milliseconds  required  for  the  internal 
pressures  to  reach  a  maximum  value.  Therefore,  the  effective  dur¬ 
ations  were  infinite  and  could  have  no  distinguishing  effect  on  animals^ 
of  different  masses;  i.  e.  ,  referring  to  Eq.  (59),  infinity  divided  by  m* 
cannot  be  distinguished  from  infinity  divided  by  m  . 

It  was  shown  above  that  the  time  parameters  of  the  blast  waves 
used  satisfied  Eq.  (70)  for  similar  animals  of  any  body  mass.  Ambient 
pressure  was  not  varied  during  the  experiments;  therefore,  according 
to  Eqs.  (69)  and  (72),  a  given  peak  blast  pressure  would  result  in  the 
same  value  of  the  index  of  biological  response  for  similar  animals 
of  any  body  mass.  Thus,  Mast  overpressure  producing  50% 
mortality  should  be  independent  of  body  mass  if  the  animal  species 
are  similar. 

The  experimental  data  from  references  (16)  and  (51)  which 
are  plotted  in  Fig.  6  suggest  that  the  animals  are  actually  divided 
into  two  groups,  one  composed  of  the  smaller  species  (mouse,  hamster, 
rat,  guinea  pig,  and  rabbit)  and  the  other  composed  of  the  larger  species 
(cat,  dog,  and  goat).  Statistical  tests  described  in  the  following  para¬ 
graph  were  made  to  determine  if  animal  mortality  within  each  of  the 
groups  is  significantly  dependent  on  body  mass. 

Assuming  that  the  P  values  plotted  in  Fig  6  were  determined 
with  the  same  reliability,  least-square  regres-  on  equations  of  the  form 

log  P  =  a  +  b  log  m 

bO 

wei  e  determined  for  the  small-animal  data  and  for  the  large-animal 
data.  The  slopes,  b,  were  examined  by  the  t-test  to  determine  the 
probability"  P^  of  obtaining  by  chance  absolute  values  of  b  as  large  as, 
or  greater  than,  those  observed  assuming  that  the  true  value  is  zero. 

The  regression  coefficients  a  and  b  are  listed  below  along  with  values 
cf  (1  -  P-D)*  the  probability  that  P^  is  dependent  on  animal  mass. 

a  b  (1  -  Pb) 

Small-animal  data:  1.  5061  .  005206  .  19 

Large-animal  data:  1.6060  . 0<8U4  .7 8 

Since,  as  shown  above,  animal  mortality  within  each  of  the 
groups  is  not  significantly  dependent  on  body'  mass,  the  P  data  were 
average  to  obtain  31.9  psi  for  the  small  animals  and  48.  f  psi  for 
the  larger  ones.  These  groupings  are  consistent  with  the  lung-volume 
and  lung-density  data  plotted  in  Fig.  2  which  indicate  that  mouse,  rat, 
guinea  pig,  and  rabbit  are  members  of  one  group  and  monkey,  cat, 
dog,  and  man,  members  of  another.  (No  data  were  shown  for  the 
hamster  and  ti  e  goat.) 

5.  4  Long-Duration  Blast  Waves  with  a  Time  Step 

1 8 

In  other  experiments  small  animals  were  placed  in  wide-mesh 
cages  and  exposed  side-on  at  the  end-plate  of  a  shock  tube  and  at  various 
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■  Reflected  Overpressure  Producing  50%  Mortality  When  The  Animal 


scmkiii;  equations  the  response  of  similar  animals  to  a  "iven  blast  wa.e  of  this 
type  would  be  identical  (Data  from  Kef  1<>  and  bj) 


distances  up  to  12  in.  away  from  it.  An  attempt  was  made  to  produce 
incident  and  reflected  shockwaves  of  constant  magnitude  by  using 
approximately  the  same  overpressure  in  the  compression  chamber. 

The  reflected  overpressures  measured  for  the  entire  series  of  shots 
varied  from  48  to  56  psi,  the  average  being  52.  1  psi.  The  corresponding 
average  of  the  incident  shock  overpressures  was  between  17  and  18  psi, 
spanning  the  17.2  psi  shock  that  theoretically  would  produce  a  reflected 
wave  of  52.  1  psi  for  the  conditions  of  perfect  reflection  at  normal  in¬ 
cidence  and  an  ambient  pressure  of  12  psi. 

Animal  mortality  vs.  distance  from  the  end-plate  were  reported 
for  the  mouse,  rat,  guinea  pig,  and  rabbit.  All  animals  placed  against 
the  end-plate  were  killed  but  the  mortality  showed  a  decline  as  the 
distance  from  the  end-plate  was  increased  --  the  distance  necessary 
to  afford  ^  given  degree  of  protection  being  greater  for  the  large  animals 
than  that  for  the  small  ones. 

Scaling  equation  (69)  was  satisfied  since  the  average  blast 
overpressures  (incident  and  reflected)  were  approximately  the  same 
for  all  species.  The  biological-response  numerics,  Eq.  (72),  would  be 
the  same  for  all  species  if  the  time  numerics,  described  in  Eq.  (70), 
were  the  same.  The  characteristic  times  significant  to  these  experi¬ 
ments  were,  of  course,  the  times  between  the  incident  and  reflected 
blast  waves.  Thus,  according  to  Eq.  (70),  (time  step)  /  (body  mass)^  3 
is  predicted  to  be  constant  for  ali  species  for  a  given  index  of  biological 
responsci  Eq.  (72).  (Note  that  ambient  pressure  was  not  varied  in 
these  experiments.  ) 

The  scaling  concepts  presented  above  were  applied  to  the 
experimental  data  by  plotting  in  Fig.  7,  for  each  species,  scaled  time 
step  vs.  body  mass  for  the  condition  of  50%  mortality.  The  average 
scaled  time  step  for  the  mopse,  rat,  guinea  pig,  and  rabbit,  was 
found  to  be  0.  516  rnsec/kg1  3,  applicable  to  the  experimental  conditions 
described.  The  plotted  data,  Fig.  7,  indicate  no  dependence  of  the 
scaled  time  steps  on  body  mass  and  thus  demonstrate,  within  reasonable 
limits,  the  similarity  of  these  species  of  animals. 

Ir.  the  preceding  paragraphs  the  scaling  of  a  time  parameter 
of  file  blast  wave  with  animal  body  mass  was  demonstrated;  however, 
nothing  was  said  in  regard  to  the  mechanisms  involved  in  reducing 
mortality  by  moving  the  animals  away  from  ti  a  reflecting  surface. 

To  help  explain  these  mechanisms,  numerical  solutions  were  obtained 
for  appropriate  blast  waves  using  the  model  equations  and  the  adjusted 
animal  parameters  listed  in  Table  1  scaled  to  a  1.7  -kg  rabbit,  the 
average  mass  of  the  rabbits  used  in  the  experiment  described  aoove. 
Model  solutions  are  shown  in  Fig.  8  for  the  situation  of  the  animal 
cage  being  against,  and  at  three  noted  distances  away  from,  the 
end-plate  which  situations  were  associated  with  '00,  80,  12,  and  0 
percent  mortality,  respectively .v 

*  The  time  between  the  incident  and  reflected  shock  in  Fig.  8(a) 
corresponds  to  the  travel  time  of  the  blast  wave  from  the  midpoint  of 
the  animal  cage  to  the  reflecting  surface  and  return. 
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Fig.  8  Computed  cavity,  or  thoracic,  o\  vrpressure  vs.  time  for 

1.  7 -kg  rabbits  exposed  at  various  distances  from  the  end-plate  of 
tile  shock  tube  The  incident,  and  refierr.eci  shocks  were  17.  2  and 
52.  1  psi,  respectively.  The  mortality  data  shown  are  from  if.ef  18. 


Note  for  the  four  blast  situations  that  the  peak  cavity  pressure 
is  inversely  related  to  the  minimum  cavity  pressure  occurring  after 
the  arrival  of  the  reflected  shock,  i.  e.  ,  cavity  pressure  resulting 
from  the  incident  wave  "inhibits"  the  effect  of  the  reflected  wave.  The 
computed  data  in  Fig.  8(c)  show  that  the  cavity  pressure  had  reached 
a  maximum  value  of  44  psi  at  the  time  the  reflected  shock  of  52.  1 
psi  arrived  and  that  the  subsequent  peak  cavity  pressure  was  only  58 
psi.  At  the  other  extreme,  the  data  in  Fig.  8(a)  indicate  that  the  cavit> 
pressure  was  4  psi  when  the  reflected  wave  arrived  which  resulted  in 
a  peak  cavity  pressure  of  333  psi. 

Two  reasons  are  proposed  to  account  for  the  reduced  mortality 
when  the  animal  is  moved  away  from  the  reflecting  surface.  (1)  The 
driving  pressure  (the  difference  between  the  external  and  internal 
pressure)  of  the  reflected  shock  is  reduced  because  of  the  increased 
pressure  in  the  lungs.  (2)  Tlje  increased  internal  pressure  also  re¬ 
sults  in  a  stronger  air  spring  being  active  at  the  time  of  arrival  of 
the  second  pulse  and  tends  to  inhibit  further  motion  of  the  chest  wall 
and  abdomen.  These  effects  are  maximized  --  the  animal  is  most 
protected  --  if  the  reflected  shock  arrives  at  the  time  when  the  lung 
pressure  has  readied  a  peak  value  due  to  the  incident  wave.  Compare 
Figs.  8  (b).  8(c),  8(d). 

5.  5  Dog  and  Goat  Blast  Data  Scaled  to  a  70 -kg  Mammal  and 
from  12. 0-  to  14.  7-  psi  Ambient  Pressure 

Reference  10  contains  a  summary  of  the  results  of  experiments 
in  which  dogs  and  goats  were  exposed  i  gainst  reflecting  surfaces'’1''  '''  to 
blast  waves  generated  in  shock  tubes  and  by  high  explosives.  The  average 
atmospheric  pressure  at  the  experimental  station  was  12.0  psi. 

Maximum  overpressures  and- the  corresponding  blast-wave  durations 
were  reported  for  the  condition  of  50%  mortality  within  24  hours  after 
exposure. 

Tolerance  to  similar  blast  situations,  but  at  sea-level  atmos¬ 
pheric  pressure,  was  cstimateu  for  a  70-kg  mammal  by  applying  the 
scaling  equations  presented  in  Secs.  5.  1  anc  5.  2  to  the  experimental 
data.  (These  would  also  be  reasonable  estimates  for  man  to  the  degree 
that  his  thoracic  structure  is  similar  to  that  of  the  experimental  animals. 
See  Sec.  4.2).  The  index  of  biological  response,  Eq.  (72),  remains 
constant  if  F.qs.  (69)  and  (70)  arc  satisfied.  Thus,  applicable  to  the 
data  for  both  dogs  and  goats, 

P'  =  P'  (14.  7/12.  0)  ( 734 

2  1 

where  P'  is  blast  pressure  for  an  ambient  pressure  of  14.  7  psi  if  P: 

3  l 

is  the  blast  pressure  determined  for  12.  0  psi. 


See  Sec .  3.2  . 

Experimental  evidence  of  the  effects  of  ambient  pressure  on 
the  tolerance  of  mice  to  air  blast  can  be  found  in  Ref.  53. 
Note  that  this  type  of  exposure  minimizes  the  wind  effect. 
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Scaled  from  the  data  for  the  16.  5-kg  dogs,  blast  wave  duration 
is  defined  by 


t 

-  t  ,  <70/16.  5)l/  3 

(12.  0/14.  7) 17  2 

(74) 

2 

id 

and  for  the 

22 .  2 -kg  goats 

t 

=  t  (70/22.  2) 17  3 

(12.014/7  )17  3 

(75) 

3  ig 

where  t2  is  duration  applicable  to  a  70 -kg  mammal  and  an  ambient 
pressure  of  14.  7  psi  if  t:  d  and  tig  are  durations  for  dogs  and  goats, 
respectively,  where  the  ambient  pressure  is  12.0  psi. 

The  original  and  the  scaled  data  are  plotted  in  Fig.  9.  The 
smooth  curve  is  an  eye-fit  of  the  data  scaled  for  the  70-kg  mammal 
and  to  sea  level  conditions,  The  consistancy  of  the  results  is  naturally 
worthy  of  comment  and  may  be  of  considerable  importance  in  estimating 
human  tolerance  to  blast  overpressures.  This  is  probably  true  if  it 
turns  out  that  man  responds  more  like  the  "larger"  than  the  "smaller" 
mammalian  species. 
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Maximum  Reflected  Overpressure,  psi 


Estimate  of  the  Parameters  of  Reflected  Blast  Waves  Resulting  in 
50%  Mortality  of  70  Kg  Mammals  at  14  7  psi  Ambient  Pressure 


Scaled  from  experimental  data  for  16.5  Kg 
dogs  of  12.0  psi  ambient  pressure 

Scaled  from  experimental  data  for  22.2 Kg 
goats  at  12  Opsi  ambient  pressure 

Original  data 


High  Explosives  —  [  —  Shock  Tube 


2  5  '10  20  50  lOO  200 

Duration,  msec 

t  -g  ■  9  An  example  oi  body-mass  and  ambi  i-nt-pr es  sure  sc, 
using  experimental  data  from  Ret  10.  The  smooth  curve  is 
eye  III  of  the  dog  ai.d  goat  data  scaled  lo  a  70-kg  mammal  ai 
n.  i  a 1 1 1  1j; e i ,i  jji'uasum.  Tim  following  readings  were  l a i 

from  the  curve: 

2  msec  290  psi 

3msec  1  8  8  p  s  1 

3  msec  1>0  psl 

10  msec  79  psi 

20  msec  54  ps, 

30  msec  52  psi 

•500  msec  52  psi 


6.  0  DISCUSSION 

6.  1  General 

In  general  it  is  well  to  emphasise  that  the  present  study 
represents  an  exploratory  effort  to  bring  theory  and  experiment  to¬ 
gether  with  the  objective,  of  trying  to  understand  more  fully  the  mech¬ 
anisms  that  are  responsible  for  primary  air  blast  injury.  The  attempt 
at  integrating  physical,  biophysical  and  biological  factors  suffers  not 
only  because  of  "state-of-the-art1’  limitations,  conceptual  and  technical, 
but  also  because  necessary  information  is  simply  not  yet  available 
even  though  some  of  it  could  he  obtained  making  use  of  appropriate 
experimental  j  ocedurcs. 

It  is  an  encouraging  matter  that  the  complex  thoracic  and 
abdominal  portions  of  the  body  can,  for  the  purpose  of  helping  visualise 
the  effects  of  blast  loading,  be  viewed  simply  as  a  gas -filled  box  to 
which  are  attached  an  orifice  and  two  pistons  to  simulate  flow  in  the 
airways  and  the  movement  ol  the  chest  wall  and  the  diaphragm, 
respectively.  It  is  not  surprising  that  the  somewhat  tedious  mathe¬ 
matical  and  analytical  formulations  were  made,  programmed,  and  that 
actual  model  solutions  were  computed.  What  is  more  significant, 
however,  is  that  considerable  progress  u.is  made  not  only  in  defining 
the  parameters  of  the  thud-mechanical  model  conceived,  bu  in  actually 
obtaining  values  that  are  reasonable  for  the  several  individual  parameters. 
Though  these  are  admittedly  approximations,  it  was  nonetneless  true 
that  only  relatively  small  adjustments  were  needed  to  obtain  satisfying 
agreement  between  computed  and  measured  thoracic  pressures  in 
the  instances  described  for  rabbits  exposed  to  shock  tube  generated 
blast  waves.  While  such  results  suggest  that  the  parameter  values 
have  a  much  better*than  order -of-rnagnitude  reliability  at  the  present 
time,  it  should  be  candidly  stated  that  there  is  much  yet  to  be  done  in 
adding  sophistication  to  improve  accuracy  in  the  approach  taken. 

The  model  described  might  have  been  formulated  using  one 
piston  rather  than  two.  Since  the  parameters  estimated  for  the  two 
pistons  were  such  that  they  act  as  one;  i.  c.  ,  their  time  displacement 
histories  were  identical  for  a  given  blast  experience.  Although  this  is 
unlikely  to  represent  the  "real-life"  situation,  presently  available 
information  does  not  allow  the  action  of  the  abdomen  to  be  differentiated 
from  that  of  the  chest  wall.  However  work  already  planned  with  the  model 
requires  a  two-pistoi,  system. 

6.  2  The  Geometric  Model  of  the  Thorax 

A  geometric  model  of  the  thorax  was  formulated  to  help 
estimate  realistic  rotations  between  the  gaseous  volume  of  the  lungs 
and  (1)  the  area  of  the  pleural  surface  and  (Z)  the  mass  of  the  chest  wall 
tissue  which  moves  in  the  blast-induced  implosion  process.  While  the 
right  circular  cylinder  with  a  height  three-quarters  of  its  diameter 
(see  Fig  3)  had  approximately  the  same  area-volume  relations  as  those 
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measured  tor  the  cadaver  used  for  this  study  it  is  dear  that  the  geo¬ 
metric  shape  chosen  is  far  from  ideal  for  application  to  all  mammalian 
species,  including  man.  In  fact,  work  is  now  under  way  employing  a 
truncated  cone  as  a  geometric  model,  a  much  more  realistic  shape. 
Also  the  analytical  details  arc  being  refined  to  allow  the  are;  s  of  the 
sides  and  the  base  of  the  cone  to  change  as  the  simulated  lung  volume 
changes  during  the  implosion  process.  Additional  r elinerncnts  may 
result  from  studies  which  are  planned  in  which  cross -sections  similar 
to  those  described  in  Sec.  4.  3  will  be  made  for  various  experimental 
animals. 


6.  3  The  Estimated  parameters  for  1 -kg  Animals  and 
Intraspecies  Scaling 


The  fact  that  physiological  data  were  presented  indicating 
that  mammalian  species  appeared  to  fall  into  two  groups  --  one  for 
''large"  and  the  other  lor  "small"  animals  --  deserves  considerable 
emphasis  ,  This  segregation  of  the  animals  of  interest  was  reflected  in 
the  best  estimates  for  model  parameters  presented  in  Table  1  wherein 
figures  were  given  for  "small"  and  "large"  species  all  scaled  to  the 
1-kg  animal  along  with  reference  to  the  scaling  equations  allowing  one 
to  estimate  the  parameters  for  an  animal  of  arbitrary  size.  Again  it 
is  necessary  to  emphasize  that  the  numerical  data  given  represent 
only  approximations,  but  it  is  easy  to  appreciate  the  desirability  and 
need  for  obtaining  accurate  numbers  for  all  parameters  along  with  a 
fair  assessment  of  the  normal  variability  associated  with  each.  For 
example,  one  of  the  problem  areas  in  blast  and  shock  biology,  that  o. 
estimating  human  tolerance  to  blast-induced  pr<  ssure  variation,  would 
be  better  illuminated  if  the  model  parameters  in  Table  1  were  known 
to  be  proper. 


6.  4  Review  of  Various  Estimates  of  Man's  Response  to 
Overpressure  ~ 

First,  consider  Fig..  10  showing  the  results  obtained  by 
Fisher,  Krohn  and  ZuckermanJ’  using  m.ice,  guinea  pigs,  rabbits, 
and  monkeys  exposed  in  wire  cages  and  goats  restrained  in  the  stand¬ 
ing  position  with  rope  netting,  tail-on  to  the  detonation  ol  1,  8,  and 
66  2/  3  lb.  PAG  charges  at  ground  level.  Pulse  durations  were  stated 
to  range  from  1  to  3  msec.  Fig.  11,  giving  the  empirically  extrapolated 
relationship  between  the  P^q  pressure  and  body  weight,  depicts  a  predicted 
Pr q  of  392  psi  for  the  60-  and  460  psi  for  the  80-kg  animal.  In  addition 
the  authors  cited  field- study  results  of  12  instances  of  human  exposure 
to  blast  for  which  the  associated  maximum  overpressures  were  estimated 
to  have  been  from  170  to  600  psi.  The  only  death  reported  was  for  a 
situation  for  which  the  maximum  overpressure  was  estimated  to  be  450  psi. 
Most  of  the  individuals  who  survived  showed  symptoms  of  blast  injury. 

Cf  these,  10  exposures  were  judged  te>  be  between  170  and  450  psi  and 
the  remaining  case  between  500  and  600  psi.  It  would  help  had  estimates 
been  reported  for  the  duration  of  the  blast  waves  and  for  the  weights 
and  tvnpfl  of  clia***7**^  t o inf'  tHp’n 
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WLLA'i  iuil  BETWEEN  (BODY  Yv  T);:/1  AND  50%  LETHAL 
HI. AT!  PRESSURE  EON  MICE,  GUINEA  FIGS  AND  RABBITS. 
HU  SiKAIOHT  LINE  IS  1  HE  CALCJLATLD  LINE 
OF  CLOSEST  FIT 
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RELATION  BETWEEN  LETHAL  EFFECT  OF  BLAST  AND 
BODY  WEIGHT  IN  MONKEYS  AND  GOATS.  THE  DOTTED  LINE 
REPRESENTS  THE  EMPIRICAL  RELATION  SHOWN  IN  THE  CHART  ABOVE 
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Fig.  10  Mortality  data  from  studies  by  Fisher,  Krohn, 
and  Zuckerman^  in  which  animals  were  exposed  tail-on 
to  the  detonation  of  high  explosive  charges  at  ground  level. 


Second,  attention  is  directed  to  Fig.  12,  a  summary  type 
illustration  taken  from  the  work  ul  Scharain,  ^1  indicating  on  a 
weight-of-char  ge  vs.  distance  diagram  that  tolerance  curves  lor 
physical  and  biologic  media  have  certain  similarities;  vu.  ,  the  curves 
are  nearly  parallel  with  the  isopressure  lines  lor  "larger"  charges 
at  the  "greater"  distances  ("long"-duration  waves)  and  approacli 
parallelism  with  the  isomomentiun  lines  lor  the  "smaller "charges 
at  the  "lesser"  ranges  ("short" -duration  waves) .  The  curve  Schardin 
cited  for  dogs,  drawn  from  the  data  of  Desaga,  ^  is  replotted  in 
Fig.  13  on  a  pressure-duration  diagram  to  show  just  fatal  conditions 
lor  dogs  which  were  exposed  on  the  ground  in  cither  recumbent  or 
erect  positions.  Open  cylinderical  charges  (height  equal  to  the  diameter) 
of  FP  60/40  (60  percent  ammonium  nitrate,  40  percent  hexogen) 
we-e  placed  and  detonated  on  the  ground.  It  is  significant  that  a  pressure 
of  215  -  220  psi  was  required  for  lethality  when  the  pulse  duration  was 
1.  b  msec,  but  when  the  latter  was  11.8  msec  the  just-fatal  pressure 
was  only  76  psi.  The  charge  weights,  lethal  limits,  pressure,  and 
duration  da^a  for  the  Desaga  experiments  arc  tabulated  in  Table  2  shown 
below. 


TAB  OF  2 


EXPLOSIVE  WEIGHTS,  DISTANCES, 
PRESSURES,  AND  DURATIONS  ASSOCIATED 
WITH  LETHALITY  IN  DOGS 
TABULATED  FROM  THE  DATA  OF  DESAGA3 


Explosive 

Weight 

kg 

Lethal 

Limit 

Meters 

Static  Maximum 
Pressure 
psi 

Pret  sure 
Duration 

i.'is  ec 

2  5 

4.  2  5 

216 

1.  0 

50 

5.  30 

-  * 

218 

1.  6 

?.on 

9.  7 

12  5 

A  T 

1000 

19.  0 

85 

8.  6 

1500 

22.  0 

7V 

10.  3 

2000 

25.  0 

76 

11.8 

That  Table  2  and  Fig.  13  demonstrate  a  deciease  in  the  lethal 
overpressure  by  almost  a  factor  of  3  (218/76  =  2.  6)  when  the  pulse  duration 
is  increased  by  more  than  a  factor  of  7  (11.  8/1.  6  -  7.  4)  deserves  consid¬ 
erable  emphasis.  Not  only  dc  tine  data  confirm  t  le  observations  of 
Hooker  ^  and  Cltuiiedsen^  that  both  pressure  and  duration  arc  important 
factors  in  blast  I  icrance  for  relatively  "short"  pulses,  but  they  represent 
the  first  lantitaiavc  results  defining  1  he  magnitude  of  the:  pres  sure -duration 
set. 


To  note  the  relationship  of  human  tolerance  for  "short"  duration 
pulses  to  the  animal  work  just  cited,  attention  is  directed  to  Fig.  14 
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DISTANCE  IN  METERS 


Fig.  12  A  charge-distance  diagram  modified  from  Schardin 
showing  characteristic  lines  of  destruction  for  animals  and 
inanimate  objects.  (Momentum  expressed  in  atm-mscc) 


Maximum  Overpressure,  psi 
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Fig.  13  Experimental  data  of  Desaga  showing  the  "lethal” 
limit"  for  dog  and  heifer.  The  point  indicated  for  man 
(100  psi,  6.  6  msec)  was  estimated  by  Desaga  as  a  result 
of  an  investigation  of  a  bombing  incident  during  the  Second 
World  War.  See  text. 
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Weight  of  Charge,  kg 


Fig.  14  (A)  Air-blast  mortality  data  as  summarized  by  BeifZinger. 

Fig.  14  (B)  Lethal  limit  data  from  studies  by  Desaga. 

Fig.  13. 


Also  see 


4  3 

showing  summarizing  charts  taken  from  Benzinger  and  Desaga. 

The  latter  shows  a  point  labeled  "man"  to  be  on  a  line  defining  lethal 
limits  tor  dugs  and  heifers  on  a  w e iglit -of -char ge  vs.  distance  diagram. 
The  former  indicates  similar  data  that  allows  one  to  appreciate  the 
relation  of  the  human  data  point  to  the  dog  experiments  cited  in 
Table  2. 

4  3 

The  human  data  referred  to  by  both  Benzinger  and  Desaga 
were  presented  in  detail  by  the  latter.  Exposure  of  8  individuals  in  an 
open-topped  anti-aircraft  gun  emplacement  occurred  following  the 
detonat  on  of  a  2,  000  lb.  HC  bomb  (total  weight  918  kg,  charge  weight 
550  kg).  All  except  one  of  the  exposed  men  were  injured,  two  fatally. 
The  fatal  cases  were  located  near  a  corner  where  pressure  reflection 
occurred.  Desaga  estimated  the  incident  overpressuie  to  be  57  psi, 
with  a  reflection  to  a  maximum  of  235  psi.  The  duration  of  the  pulse 
was  not  stated,  but  was  likely  to  have  been  about  7  msec.  Also,  it 
is  not  possible  now  to  say  how  much  tune  might  have  elapsed  between 
the  occurrence  of  the  incident  and  reflected  spike  at  the  location  of 
the  fo.t?'ly  injured  men,  but  there  might  well  have  been  a  delay  of 
more  than  1  or  2  msec,  a  fact  that  could  have  increased  tolerance 
significantly.  Finally,  one  must  point  out  that  one  of  the  fatally  injured 
individuals  aisc  suffered  a  fractured  skull,  a  fact  that  indicates  decel- 
erative  impact  as  well  as  primary  blast  might  have  beer,  a  factor  in 
lethality.  Even  so,  the  condition  of  the  lungs  described  was  typical 
of  primary  blast  injury  and  no  doubt  the  latter  contributed  mostly  to 
death  which  occurred  about  45  minutes  after  exposure  fr  r  each  case, 


Making  use  of  the  human  data  described  above,  Desaga 
estimated  that  the  "lethal  limit"  for  man  was  6,  8  atmospheres  overpressure 
(100  psi  for  an  ambient  pressure  of  14.  7  psi)  for  a  charge  of  550 
kg.  A  duration  of  6.  6  msec  was  estimated  for  this  situation  by  graphical 
interpolation  of  Desaga's  blast  data  presented  in  Table  2. 


Third,  Richmond  et  al.  .  using  data  obtained  at  Albuquerque 
(12  psi  ambient  pressure)  partly  shown  in  the  pressure  duration  diagram 
in  Fig.  15  tor  over  3,000  animals  distributed  among  6  mammalian 
species,  extrapolated  the  findings  using  both  "small"  and  "large"  animals 
to  predict  1,  50,  and  99  percent  Jcthai  conditions  for  the  70  kg  mammal 
as  noted  in  Fig.  16.  Animals  were  either  exposed  side-on  against  the 
enu-plate  of  a  shock  tube  or  lying  on  the  ground  with  charges  detonated 
vertically  above.  In  each  case,  the  animal  was  exposed  to  the  "fast"- 
rising  incident  pulse  followed  almost  immediately  by  the  reflected  pressure. 


The  approach  used  by  Richmond  et  al.  to  extrapolate  the  animal 
data  to  the  70  kg  mammal  is  illustrated  in  Fig.  17  for  "last" -ris  ing, 
400-msec  waves.  Note  that  this  figure  also  includes  data  for  the  hamster 
and  the  car  which  became  available  after  the  original  estimates  were 
made.  The  additional  data,  however,  did. not  significantly  change  the  p 
estimate  for  tr.e  70  kg  mammal,  52  psi,  which  figure  scaled  to  sea-level 
ambient  pressure  becomes  52  (14.7/12)  psi  =  04  psi.  The  latter  is  close 
to  the  62  psi  value  noted  in  Fig.  9  of  the  present  study,  a  number  obtained 
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Ovorprcssurf  For  50  P«f  Cent  LelFolitj  As  A  Function  Of  Durotion 
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16  Results  of  body-mass  extrapolation  of  the  data  shown  in  Fig.  15 
mammal,  fi  mbit'll l  pressure:  1=1  psi. 


:  Reflected  Overpressure  Producing  50%  Mortality  When  The  Animal 
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employing  the  scaling  equations  and  the  empirical  data  of  Richmond  ^ 

for  dogs  and  goats.  i 


Two  horizontal  lines  were  drawn  in  Fig.  17,  one  through 
the  "small"-  and  the  other  through  the  "large"- animal  data,  to  raise 
the  question  as  to  whether  or  not  the  animals  are  actually  divided  into 
two  distinct  groups,  one  of  which  may  include  man  (see  Sec.  5.3).  The 
alternate  assumption  is  that  all  the  experimental  animals  are  members 
of  the  same  group  and  that  tolerance  to  400-msec  blast  waves  is  a 
function  of  body  mass  as  illustrated  by  the  dashed  regression  line  in 
Fig..  17.  Although  the  scaling  or  extrapolation  procedure  used  may  not 
be  too  important  for  the  "longer" -duration  pulses,  it  is  for  the  "shorter"- 
duration  blast  waves.  Table  3  is  included  to  compare  the  results  of 
the  scaling  apprbach  described  in  Sec.  5.  5  using  only  the  goat  and  jq 
dog  data  with  the  results  of  the  extrapolation  made  by  Richmond  et  al. 
employing  both  "large"-  and  "small" -animal  data.  This  table  also 
includes  the  estimates  of  man's  tolerance  to  air  blast  made  by  Fisher, 
Krohn,  and  Zuckerman^  and  Desaga.  3 

The  information  presented  in  Table  3  is  shown  graphically 
in  Fig.  18.  Since  the  animal  work  of  Fisher  et  al.  applied  to  1  to 
3  msec  pulses  and  their  estimates  were  made  for  60  to  80  kg  animals, 
the  extrapolated  data  are  shown  in  Fig.  18  as  a  rectangle  (dashed  lines) 
whose  sides  were  placed  at  the  appropriate  limits.  It  is  apparent  that 
the  estimated  tolerance  curves  of  the  present  study  would  pass  through 
this  rectangle  if  they  were  extended  to  shorter  durations.  Desaga's 
estimate  of  the  lethal  limit  for  man,  shown  as  a  point  in  Fig.  18,  is 
located  --  perhaps  fortuitously  --  within  about  2  psi  of  the  estimated 
tolerance  curve  of  the  current  study  for  14.  7  psi  ambient  pressure.  For 
blast  waves  with  durations  greater  than  100  msec  which  are  particularly^ 
apropos  of  nuclear  weapons,  the  extrapolated  curves  of  Richmond  et  al. 
and  the  scaled  curves  of  the  present  investigation  agree  within  10  percent. 
For  durations  of  3  msec,  however,  the  former  P  estimates  are  2-3 
times  as  large  as  those  of  the  latter. 

6.  5  Other  Applications  of  the  Model 
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In  addition  to  the  use  of  the  model  described  in  the  present 
study  to  help  formulate  human  tolerance  criteria  from  data  obtained 
with  various  species  of  animals,  other  applications  are  of  interest.  For 
instance,  the  computation  of  intrathoraci.c  pressure-time  variations  in 
various  animals  exposed  in  comparable  blast  environments  not  only 
helps  one  understand  the  contributions  of  critical  body  parts  involved 
in  the  implosion  process,  but  can  also  contribute  to  the  study  of  the 
biophysical  mechanism  etiologically  active  in  causing  the  highly  hazard¬ 
ous  blast  injury  to  the  lungs. 

The  model  allows,  for  a  given  animal,  a  comparison  of  the 
effects  of  one  blast  situation  with  that  of  another  and,  for  a  given  blast 
situation,  a  comparison  of  the  response  of  one  animal  with  that  of  another. 
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TABLE  i 


A  SUMMARY  OF  VARIOUS  ESTIMATE  OF  MAN'S  TOLERANCE  TO  AIR  BLAST 

Tabulated  overpressures  in  psi  correspond  to  50  percent  mentality  except 
where  noted.  (P  :  ambient  pressure) 


- r 

P 

o 

14.7  ps. 

P 

O 

12  psi 

Duration, 

ms  ec 

Present 

Study 

Ref.  10 

Ref.  54  Kef.  3 

Present 

Study 

Ref.  10 

1  -  3 

* 

*  340 

>1000* 

392-409** 

- 

2 

c 

- 

- 

- 

3 

188 

528 

1  70 

431 

5 

120 

227 

106 

185 

6.  6 

98 

16b 

100  1 

86 

135 

10 

79 

120 

68 

98 

20 

64 

87 

1 

64 

72 

30 

62 

78 

51 

64 

60 

Cd 

71 

51 

58 

400 

62 

04 

51 

i 

52 

*  Extrapolated  to  1.??  msec. 

vv  Lower  and  higher  overpress  *  es  apply  to  oO-  and  80-hg  animal, 
respectively. 

"Lethal  limit"  values  of  overpressure  and  duration. 
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Maximum  Reflected  Ovepressure,  psi 


Duration,  msec 

Fig.  18  A  comparison  of  various  estimates  of  man's  tolerance  to 
air  blast.  See  Table  3. 
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L  'u‘  "lunr,u  "iu«h:l,  along  wit.il  the  sealing  relations  presented 
m  ihis  report,  could  be  used  in  related  studies  such  as  tin*  biological 
citccts  of  explosive  decompression,  high  noise  fields  (oscillating  ambient 
pressures),  and  a  whole  variety  of  other  disturbed  or  atypical  blast- 
induced  wave  forms,  In  some  of  these  areas,  unlike  the  case  for  the 
very  rapid  changes  in  blast  pressures  associated  with  typical  or  near- 
typical  wave  forms,  flow  in  the  airways  may  be  an  important  mechanism 
ot  compensation.  It  certainly  can  be  in  explosive  decompression, 
and  thus  a  more  accurate  air  flow  expression  than  the  one  reported 
here  woe  be  desirable. 


In  closing  the  discussion,  it  is  well  to  emphasize  tha„  the 
present  study  helped  elucidate  the  very  important  part  played  by  the 
non-linear  air  spring  which  is  estimated,  for  small  species  exposed 
at  sea-level  ambient  pressures,  to  be  mere  than  500  times  as  "stiff" 
as  the  mechanical  spring  (Kj  and  K^)  discussed  in  Sec.  4.4  and 
alluded  to  in  Eq.  (65)  and  in  Table  1.  This  significant  fact  should  not 
be*  lost  sight  of  and  it  needs  hr  explained  that  the  air -spring  effects 
enter  the  model  equations  somewhat  indirectly  through  the  pressure- 
volume  relationship  as  discussed  in  connection  with  Eqs.  (38)  -  (41). 
Thus,  the  response  time  of  the  thoracic-abdominal  system  to  blast 
loading  is  a  function  of  the  magnitude  of  the  pressure;  i.  e.  ,  the 
thoracic  pressure  will  rise  to  a  maximum  earlier  for  a  high  than  it 
will  for  a  low  overpressure.  That  high  overpressures  are  necessary 
at  hig-her  ambient  pressures  to  achieve  the  same  biological  response 
can  be  related  to  the  fact  that  the  air  spring  is  stronger  in  the  latter 
case.  Also,  the  stiff  air  spring  no  doubt  is  a  very  important  factor 
in  the  "adaptation"  an  animal  makes  to  a  non-lethal  incident  pulse  which 
completely  protects  him  against  the  associated  and  otherwise  lethal 
reflected  pulse  provided  the  two  increases  in  external  pressure  are 
separated  enough  in  time  to  allow  the  intrathoracic  pressure  to  rise 
significantly,  or  to  a  maximum,  in  response  to  the  initial  pulse.  In 
such  circumstances,  the  internal  pressures  are  high  and  thus  the 
animal  is  a  "stiffer"  system  which  can  better  resist  the  effect  of  the 
second  pulse.  This  probably  means  that  if  the  stiff  air  spring  factor 
were  not  at  work,  the  second  component  of  a  stepwise  increase  in 
air  pressure  would  be  significantly  more  effective  than  it  is. 
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SUMMARY 


7.  0 


The;  highly  lethal  effects  of  exposure  to  "fast1'  -  rising  blast- 
induced  overpressure,  the  extraordinary  sensitivity  of  the  mammal 
to  small  changes  in  the  average  rale  of  pressure  rise,  and  the  recent 
speculations  covering  the  etiologic  and  "protective"  significance  of 
the  internal  pressure  rise  as  a  consequence  of  the  violent  implosion 
of  the  body  wall  hy  blast  waves  were  cited  as  reasons  foi  exploring 
the  utility  oi  a  simplified  model  of  the  mammalian  thoraco-abdominai 
system  to  guide  thinking  about  the  implosion  process  and  its  etiologic 
relations  to  air  blast  injury. 

The  mathematical  model,  conceived  to  simulate  mammalian 
response,  incorporated  (a)  the  equations  of  motions  of  "chest-wall" 
and  "abdominal"  pistons  allowed  to  move  into  a  confined  volume  represent¬ 
ing  the  gas  contained  in  the  lung  and  airways  and  (b)  an  expression 
describing  gas  flow  through  an  orifice  serving  as  a  communication 
between  the  confined  "lung"  gas  and  that  in  the  external  environment. 

The  loilowing  parameters  of  the  model  w erL  defined: 

a.  The  average  gaseous  volume  of  the  lungs. 

b.  The  effective  areas,  masses,  spring  constants,  and 
damping  factors  ol  the  pistons  representing  the  abdomen 
and  the  chest  wall. 

c.  The  effective  area  of  the  orifice  representing  the  airways. 

The  assumptions  made  in  formulating  the  model  equations 
(which  wore  integrated  by  numerical  methods  making  use  of  a  digital 
computer)  were: 

a.  1  he  pres  sure -volume  changes  in  the  cavity  representing 
the  lungs  were  polytropic  with  an  exponent  of  1.2. 

b.  1  hough  the  gas  confined  in  abdominal  organs  obviou=’*v 
influenced  the  implosion  process,  the  effect  was  insignificant. 

c.  The  complex  movement  of  gas  to  and  from  and  within 
the  airways  of  the  lung  could  be  adequately  simulated 
with  a  simple  orifice  with  an  appropriately  chosen  coeffi¬ 
cient. 

d.  The  temperature  and  pressure  of  the  gas  in  the  lungs  were 
uniform  over  the  entire  volume  at  any  g»ven  instant. 

Scaling  equations  were  developed  for  the  animal  parameters 
based  on  the  principle  of  (a)  similarity  of  shape  and  (b)  equivalence  of  the 
tissue-mass  and  tissue-elasticity  distributions.  By  dimensional  analysis. 
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the  relation  between  body  muss  ul  similar  animal!,  and  blast-wave 
parameters  were  formulated  ler  the  condition  of  equal  biological  response, 
as  indicated  by  any  nondimensional  index,  such  as  percent  mortality. 

The  parameters  of  the  model  were  evaluated  through  a 
variety'  ol  approaches,  some  ol  which  are  noted  below. 

a.  Evidence  cited  Iroin  the  literature  indicated  that  the 

ratio  ol  the  average  lung  volume  to  average  body  mass  was 
about  9  cm  ' /kg  for  "smaller"  animals  (mice,  rats,  guinea 
pigs,  and  rabbits),  but  was  significantly  higher,  near 
30  cm0 /kg,  lor  "larger"  animals  (monkeys,  cats,  dogs,  man), 

b.  The  average  lung  densities  lor  the  "smaller"  species  were 
reported  to  be  higher  than  lor  the  "larger"  animals,  0.  367 
compared  with  0.  194  g/cm3. 

c.  A  geometric  model  of  the  thorax  was  devised  using  cross- 
sections  of  a  human  cadaver  reported  in  the  literature. 

Assuming  similarity  of  shape,  scaling  equations  were 
developed  relating  areas  and  masses  of  the  chest-wall  and 
abdominal  pistons  to  lung  volume  and  density. 

d.  Published  velocity  resonant  frequency  data  lor  cats  and 
dogs  were  used  to  estimate  the  total  spring  constant,  which 
was  apportioned  to  each  of  the  pistons  on  the  basis  ol  mass. 

e.  The  effective  orifice  area  was  scaled  from  data  reported 
for  rapid  decompression  studies  in  man. 

f.  The  damping  factor  was  estimated  by  analy/.ing  two  records 
of  intrathoracic  pressure  for  1. 8-kg  rabbits  exposed  to 
blast  waves  generated  by  shock  tubes. 

All  the  animal  parameters  were  tested  against  the  experi¬ 
mental  records  mentioned  above  (f)  by  solving  the  model  equations  for 
the  measured  shock-tube  waves  using  parameters  scaled  to  a  I.  8-kg 
rabbit.  Adjustments  were  made  in  most  of  the  parameters  in  order  to 
achieve  better  agreement  between  the  computed  and  measured  thoracic 
pressures.  Graphical  comparisons  were  made  between  the  experimental 
and  computed  pressure  records  and  the  adjusted  parameters  were  listed 
along  with  those  estimated. 

Certain  selected  applications  of  the  model  and  scaling  equations 
to  problems  in  blast  and  shock  biology  were  worked  out.  Among  them, 
along  with  other  matters  of  interest,  were  the  following: 

a.  Application  of  the  scaling  equations  indicated  that  the 

mortality  of  similar  animals  would  be  the  same,  regardless 
of  body  mass,  when  exposed  under  similar  conditions  to  a 
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1  'lung  "  -  du  r -i  1 1  un  shoe  V  Wiivi:  that  does  not  drray  appl  eciabiy 
durmK  the  short  tune  necessary  lor  1 1  e  lung  pressures  to 
reach  the  peak  values.  A  test  ol  tins  theory,  using  published 
data,  indicated  that  the  mouse,  hamster,  rat,  guinea  pig,  and 
rabbit  are  approximately  similar,  hut  tli.il  the  cat,  dog,  and 
goat  have  a  signilicantiy  higher  blast  tolerance. 

c.  Published  data  were  used  to  show  (according  to  the  scaling 
equations  for  similar  animals)  that  the  mouse,  rat,  guinea 
pig.  and  rabbit  have  approximately  equal  sen:  ltivity  to  time- 
delay's  between  the  incident  and  the  reflected  shockwaves  if 
the  tunes  are  made  proportional  to  the  cube  root  of  the  body 
mass.  The  time-step  «Mart  was  also  investigated  by  solving 
the  model,  with  parameters  for  a  1.  7-ky  rabbit,  tor  blast 
waves  of  increasing  time  steps  similar  to  those  measured 
in  the  shock  tube.  It  was  found  that  internal  pressure 
resulting  from  the  incident  wave  tended  to  inhabit  the  effects 
of  flu*  reflected  wa'"’. 

d  By  scaling  procedure.',,  estimates  were  made  oi  the  para¬ 
meters  of  reflected,  "fast" -rising,  blast  waves  which  would 
result  in  50  percent  mortality  for  70-kg  mammals  when  ex¬ 
posed  at  14.  7-psi  ambient  pressure.  The  published  data  used 
to  make  these  estimates  were  results  of  shock-tube  and 
high -explosive  studies  at  12-psi  ambient  pressure  with  22.  2-kg 
goats  and  16.  5-kg  dogs  exposed  against  reflecting  surfaces. 

In  discussing  the  results,  mention  was  made  of  the  fact  that 
information  currently  available  does  notallow  the  responses  of  the  abdomen 
to  be  distinguished  from  those  of  the  chest  wall.  Work  likely  to  be  appropriate 
and  relevant  to  improvement  of  the  geometric  model  of  the  thorax  (  a  truncated 
cone  rather  than  a  cylinder)  was  noted  and  the  tentative  and  approximate 
nature  of  the  model  parameters  evaluated  for  the  idealized  1-kg  "small" 
and  "large  animal  was  emphasized. 

Three  approaches  to  assessing  man's  response  to  overpressure 
were  reviewed  and  compared  with  that  presented  in  the  present  aluuv.  Two 
of  the  estimates  for  tolerance  to  short-duration  waves  were  shewn  to  be 
cons istent  with  those  presented  herein.  The  third  approach  resulted  in 
estimated  human  tolerances  which  were  approximately  the  same  as  those 
proposed  in  this  report  for  blast  waves  of  "long"  duration  but  were 
significantly  different  for  pulses  of  "short"  duration. 
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1 3  ABSTRACT 

A  mathematical  model  was  described  which  was  developed  to  compute  socie  of  the 
fluid-mechanical  responses  of  the  thoraco-abdominal  system  subjected  to  rapid 
changes  in  environmental  pressure.  The  approach  —  helpful  in  understanding  many 
of  the  primary  effects  of  air  Mast  on  animals,  but  applicable  to  related  problems 
involving  rapid  changes  in  environmental  pressure  .s  well  --  incorporated  an  air 
cavity  representing  the  gaseous  volume  of  the  lungs  (although  it  is  recognized  that 
gas  In  the  organs  of  the  abdomen  may  Influence  the  response  of  the  system),  two 
movable  pistons  and  an  orifice  through  which  gas  might:  pass  In  either  direction. 

One  of  the  pistons  represented  the  chest  wall  and  the  other  that  portion  of  the 
abdomen  which  moves  with  the  diaphragm  to  change  the  lung  volume.  Each  piston  was 
"assigned"  an  effective  mass  and  area,  a  spring  constant,  and  a  damping  factor. 

The  orifice  was  taken  to  "incorporate"  the  characteristics  of  the  many  airways  of 
the  respiratory  system.  Parameters  relating  the  animal  to  the  model  were  estimated, 
tested  and  then  adjusted  as  required  by  comparing  model  results  with  experimental 
records  of  thoracic  pressures  recorded  for  rabbits  exposed  to  blast  waves  in  shock 
tubes . 

Equations  were  derived  to  Beale  parameters  applicable  to  a  given  animal  to 
those  for  similar  creatures  of  arbitrary  mass.  By  dimensional  analysis  other 
equations  were  developed  to  relate,  for  a  given  biological  response,  the  body  mass 
of  similar  animals  to  blast  wave  parameters.  Numerical  solutions  of  the  model  were 
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presented  to  help  explain  the  mechanisms  involved  when  animals  were  "loaded" 
with  typical  wave  forms  or  with  pulses  increasing  to  a  maximum  in  a  stepwise 
manner,  a  contingency  associated  with  a  quite  significant  increase  in  mammalian 
tolerance  to  over-pressure.  Differences  in  response  to  "short-"  and  "long"- 
du rat ion  blast  waves  were  noted.  Applications  of  the  scaling  concepts  were 
exemplified  in  several  ways  ma.’:lng  use  of  the  published  data  in  blast  biology. 
In  one  instance,  the  blast  tolerance  of  a  70-kg  mammal  was  estimated  for  sea- 
level  ambient  pressure  making  use  of  experimental  data  for  dogs  and  goats 
obtained  at  Albuquerque  altitude  (ambient  pressure  of  12  psi).  That  the  latter 
might  have  significant  implications  in  assessing  human  response  to  blast-pro¬ 
duced  overpressures  was  discussed  along  with  several  other  relevant  matters. 


